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Chapter  1 


Introduction 

1.1  Background 

It  is  of  interest  to  be  able  to  predict  and  study  the  radar  cross  section  (RCS)  of  targets  on  a 
rough  sea  surface  using  computational  electromagnetics  (CEM)  techniques.  Figure  1.1  shows 
a  typical  situation  where  a  ship  is  illuminated  by  an  incident  electromagnetic  (EM)  plane 
wave,  and  scatters  energy  in  all  directions.  It  is  known  that  the  prevailing  sea  state  (surface 
roughness)  can  greatly  affect  the  radar  scattering  of  a  ship.  The  time-changing  sea  surface 
may  increase  or  decrease  the  amount  of  energy  backscattered  to  the  receiver  compared  to  a 
calm  surface,  so  the  RCS  likewise  changes  with  time.  To  predict  the  effect  of  the  sea  surface 
on  the  RCS,  it  is  necessary  to  study  how  the  incident  field  scatters  from  the  rough  sea  surface 
before  hitting  the  ship,  and  how  the  subsequent  scattering  from  the  ship  interacts  with  the 
sea  surface  before  propagating  back  to  the  receiver. 


Figure  1.1:  Scattering  by  a  ship  on  a  rough  sea  surface  illuminated  by  an  incident  EM  plane 
wave. 
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This  report  describes  the  research  performed  during  the  first  year  of  Grant  No.  N00014- 
98-1-0243  from  the  Office  of  Naval  Research.  The  particular  focus  is  to  develop  algorithms 
wherein  existing  CEM  tools  for  predicting  the  RCS  of  a  target  can  be  modified  to  account 
for  the  presence  of  the  rough  sea  surface,  without  having  to  include  the  infinitely  large 
rough  surface  in  the  computational  domain.  The  first  year’s  study  is  devoted  to  a  thorough 
analysis  in  two  dimensions  (2D)  of  the  phenomenology  associated  with  the  target/rough 
surface  scattering  problem.  Basic  algorithms  are  developed  and  tested  to  provide  proofs-of- 
concept  before  proceeding  to  more  complex  3D  problems.  It  is  noted  that  the  backscatter 
from  the  sea  surface  itself  (i.e.,  “sea  clutter”)  is  not  of  interest  here,  but  has  been  studied 
extensively  [1,  2].  It  is  also  noted  that  plane  wave  incidence  is  of  primary  interest  here 
because  the  definition  of  RCS  assumes  that  the  source  is  at  infinity.  For  typical  surface-to- 
surface  or  air-to-surface  radar  search  scenarios,  this  assumption  implies  that  the  downward 
propagating  incident  field  is  locally  planar  in  the  vicinity  of  the  target  (but  the  incident  field 
reflected  from  the  sea  surface  is  not  assumed  to  be  locally  planar).  The  propagation  of  the 
EM  field  from  the  source  to  the  target,  as  well  as  the  source  and  receiver  gain,  is  not  included 
in  the  plane  wave  RCS  but  is  included  elsewhere  in  the  standard  radar  equation  [2],  The 
power  received  by  a  radar  antenna  with  effective  aperture  Ae  due  to  scattering  from  a  target 
with  RCS  a  is  given  by 

Pr  =  Pt  G  Ae  S2  a  (1.1) 

where  Pt  and  G  are  the  power  and  gain  of  the  transmitting  antenna,  and  S  is  the  one-way 
propagation  loss  factor.  For  free  space  S  is  1/4ttR2,  where  R  is  the  range  to  the  target.  For 
propagation  over  the  sea  surface,  S  is  a  much  more  complicated  function  of  distance  and  is 
not  the  subject  of  this  report  (although  the  GFB  method  of  Chapter  2  can  compute  this 
function  for  limited  propagation  distances). 
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Incident  Field 

Figure  1.2:  Ray-optical  characterization  of  the  field  illuminating  the  ship. 

1.2  Technical  Approach 

Since  most  existing  CEM  tools  for  RCS  prediction  rely  on  a  ray-optical  incident  field,  it  is 
of  primary  importance  to  be  able  to  characterize  the  sea  scattered  field  which  illuminates 
the  ship  ray-optically,  as  shown  in  Figure  1.2.  This  is  especially  important  for  RCS  tools 
which  use  high-frequency  asymptotic  methods,  such  as  physical  optics  (PO),  the  physical 
theory  of  diffraction  (PTD),  and  the  uniform  theory  of  diffraction  (UTD)  [3].  These  methods 
must  have  a  well-defined  ray-optical  incident  field.  Furthermore,  most  ship  geometries  are 
very  large  in  terms  of  the  wavelength  of  the  EM  field  at  radar  frequencies,  so  high-frequency 
methods  must  usually  be  used  instead  of  numerical  methods  which  are  better  suited  for 
smaller  problems. 

The  ray-optical  characterization  investigated  here  uses  a  discrete  plane  wave  expansion  for 
the  sea  scattered  incident  field,  as  shown  in  Figure  1.3.  Once  the  incident  field  is  characterized 
ray-optically,  the  scattering  from  the  ship  is  found  approximately  by  replacing  the  rough  sea 
surface  with  a  locally  planar  surface,  as  shown  in  Figure  1.4.  This  allows  the  half-space 
Green’s  function  to  be  used  to  compute  the  scattered  field,  or  equivalently,  image  theory 
could  be  used  [4],  (Although  image  theory  is  not  recommended  because  including  the  image 
of  the  ship  doubles  the  size  of  the  problem.)  The  ray-optical  expansion  of  the  incident  field 
combined  with  the  local  planar  approximation  for  the  sea  surface  effectively  decouples  the 
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target  scattering  problem  from  the  sea  surface. 

The  ship  scattered  field  is  found  only  over  the  surface  of  the  ship.  With  this  information, 
the  reciprocity  theorem  [4]  is  used  to  compute  the  backscattered  field  at  the  source  location. 
The  reciprocity  theorem  is  an  exact  formulation  which  gives  the  target  scattered  field  at  the 
source  location  in  terms  of  an  integral  over  the  surface  of  the  target.  The  integral  is  the 
reaction  of  the  total  incident  field  (including  the  incident  field  scattered  from  the  sea  surface) 
with  the  target  scattered  field.  The  only  approximation  is  the  use  of  a  locally  planar  surface 
in  the  vicinity  of  the  target  to  find  the  scattered  field.  The  use  of  reciprocity  eliminates  the 
need  to  track  the  scattered  fields  back  to  the  source. 

Since  an  approximation  is  used  to  find  the  target  scattered  fields,  the  decoupled  ray- 
optical  RCS  results  will  not  be  in  perfect  agreement  with  the  reference  solution  which  con¬ 
siders  the  sea  surface  and  target  as  a  whole.  However,  the  goal  of  this  research  is  not 
necessarily  to  predict  the  exact  RCS  (which  changes  with  time  as  the  sea  surface  evolves), 
but  to  predict  the  peak  levels  and  variations  in  the  RCS  for  a  given  sea  state. 

1.3  Organization  of  Report 

Chapter  2  presents  an  efficient  numerical  approach  termed  the  Generalized  Forward-Backward 
method  for  computing  the  scattering  from  a  2D  target  on  a  rough  sea  surface.  This  method 
is  used  throughout  the  report  to  generate  reference  solutions  for  the  ray-optical  algorithms, 
and  to  study  the  scattering  properties  of  different  targets  under  various  conditions. 

Chapter  3  studies  the  plane  wave  spectral  properties  of  the  sea  scattered  field  with  the 
target  absent.  The  coherent  (time-independent)  and  incoherent  (time-varying)  components 
are  clearly  delineated.  A  discrete  plane  wave  expansion  for  the  field  which  illuminates  the 
target  is  then  proposed.  Chapter  4  then  uses  this  expansion  as  the  excitation  for  finding  the 
scattering  from  the  target  in  the  presence  of  a  locally  flat  surface,  effectively  decoupling  the 
target  scattering  problem  from  the  sea  surface.  The  reciprocity  integral  is  formulated  and 
used  to  compute  the  RCS  (backscattered  field)  by  integrating  only  over  the  target  surface. 

Chapter  5  investigates  the  behavior  of  the  RCS  as  a  function  of  the  time-changing  sea 
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surface,  and  as  a  function  of  the  roll  angle  of  the  target.  Different  target  geometries  are 
studied,  including  a  low-RCS  shape.  These  results  give  very  good  insight  into  the  scattering 
phenomenology,  and  indicate  the  usefulness  of  the  ray-optical  techniques  for  predicting  peak 
levels  and  the  pattern  variation  of  the  RCS.  Finally,  conclusions  and  future  directions  of 
research  are  discussed  in  Chapter  6. 
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Chapter  2 


The  Generalized  Forward- Backward 
Method 

2.1  Introduction 

In  this  chapter  a  generalization  of  the  Forward-Backward  (FB)  method  is  presented.  This 
generalization  allows  us  to  study  the  scattering  from  composite  surfaces  that  can  include 
one  or  more  large  arbitrarily  shaped  obstacles  (like  a  ship  or  a  large  rogue  breaking  wave) 
on  the  ocean  surface,  as  shown  in  Figure  2.1.  The  conventional  FB  method  is  not  expected 
to  exhibit  convergent  behavior  for  such  problems.  The  new  approach,  called  the  Generalized 
Forward- Backward  (GFB)  method,  is  based  on  a  combination  of  the  conventional  FB  method 
with  the  Method  of  Moments  (MoM),  where  the  MoM-  is  only  applied  to  the  region  close  to 
the  obstacle.  The  solution  is  found  through  an  iterative  procedure  based  on  the  same  general 
concepts  as  the  FB  method,  but  with  some  significant  differences.  The  computational  cost  of 
the  GFB  method  is  similar  to  the  FB  solution.  It  only  includes  an  additional  cost  associated 
with  the  direct  MoM  solution  of  a  small  region  containing  the  obstacle  and  nearby  sea 
surface.  The  GFB  method  is  very  useful  for  predicting  and  studying  the  EM  scattering  from 
targets  in  the  presence  of  a  rough  surface.  It  also  used  here  to  provide  reference  solutions  to 
validate  the  more  approximate  ray-optical  approaches  developed  in  later  chapters. 

This  chapter  is  organized  as  follows:  the  conventional  FB  formulation  for  the  solution  of 
the  electric  field  integral  equation  (EFIE)  is  presented  in  Section  2.2.  Section  2.3  develops  the 
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Figure  2.1:  A  scattering  obstacle  on  an  ocean-like  surface  illuminated  by  a  low-angle  incident 
field. 

generalization  of  the  previous  method  which  leads  to  the  new  GFB  iterative  solution.  Some 
numerical  results  and  a  summary  and  conclusions  are  also  included.  A  Spectral  Acceleration 
of  the  Generalized  Forward-Backward  method  is  presented  in  Section  2.4.  In  the  rest  of  the 
chapter,  all  the  fields  and  the  currents  will  be  assumed  to  have  a  time-harmonic  dependence 
of  the  form  eJwt  which  will  be  suppressed  from  the  field  expressions.  The  radian  frequency 
is  u,  and  n  and  k  are  the  permeability  and  wavenumber,  respectively,  of  the  medium  above 
the  rough  surface  (generally  assumed  to  be  free  space). 

2.2  The  Conventional  Forward- Backward  Method 

The  magnetic  field  integral  equation  (MFIE)  for  a  two-dimensional  scattering  problem  deal¬ 
ing  with  a  one-dimensional  rough  surface  and  its  solution  by  the  FB  method  has  been 
presented  in  [5].  Here  we  are  concerned  with  the  solution  of  the  EFIE  in  a  similar  way. 
The  application  of  the  FB  method  combined  with  the  EFIE  has  already  been  done  in  [6], 
where  a  new  approach  based  on  the  spectral  representation  of  the  Green’s  function  acceler¬ 
ates  the  computation  of  the  forward  and  backward  matrix-vector  products.  In  this  section 
the  formulation  of  the  FB  method  is  briefly  reviewed,  starting  with  the  EFIE  followed  by  a 
discretization  process  using  the  MoM  to  establish  the  corresponding  matrix  equations  which 
will  be  solved  in  an  iterative  way. 
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Figure  2.2:  (a)  Ocean-like  surface  illuminated  by  a  TM  ( Ey )  polarized  near-grazing  incident 
field,  (b)  Forward  and  backward  regions  for  the  n-th  matching  point. 


Consider  a  one-dimensional  ocean-like  surface  depicted  in  Figure  2.2(a).  The  sea  water 
is  modeled  as  a  perfect  electrical  conductor  (PEC),  so  there  is  no  penetration  into  the  sea 
surface  S.  (The  method  may  also  be  applied  to  impedance  surfaces  without  loss  of  generality.) 
The  horizontal  and  vertical  coordinates  of  this  parametric  surface  are  respectively  x(£)  and 
z(£),  both  functions  of  the  path  length  l  on  the  surface. 

For  such  a  surface,  illuminated  by  a  TM  ( Ey )  polarized  incident  field,  the  EFIE  can  be 
obtained  by  imposing  the  tangential  electric  field  boundary  condition  on  the  PEC  surface  S : 

~  jsJ(p')  (k  |  p  -  p'|)  df  =  - E‘(p )  (2.1) 
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where  El{p)  is  the  incident  field  at  a  point  p  €  S,  J(p')  is  the  induced  current  and 
( k  |p  —  p'|)  is  the  second  kind  Hankel  function  with  order  zero  [7].  The  current  J(p')  ra¬ 
diates  in  free  space  and  generates  the  scattered  fields.  This  equation  is  to  be  solved  through 
discretization  by  using  the  MoM.  Although  the  surface  S  extends  from  minus  infinity  to 
infinity  in  x,  the  incident  field  is  considered  to  be  tapered  so  that  the  illuminated  rough 
surface,  and  thus  the  integration  in  (2.1),  can  be  confined  to  a  finite  region  of  length  L. 

For  our  purposes,  we  will  apply  the  most  simple  formulation  of  the  MoM  using  a  set  of  N 
pulse  basis  functions  and  point  matching  weighting  at  the  center  of  each  current  element  [7]. 
Typically,  about  10  pulse  basis  functions  per  wavelength  are  used.  After  the  discretization 
process,  equation  (2.1)  is  transformed  into  a  matrix  equation: 

1-7  =  7  (2.2) 


The  elements  of  the  impedance  matrix  Z  are  given  by: 


Zrnn 


■<fH0 


(2) 


(k  \Pm~Pn\)Ai 


m  =  n 
m^n 


(2.3) 


where  7  is  the  Euler  constant  1.871,  Al  is  the  width  of  the  pulse  basis  and  pn  is  the  position 
vector  of  the  nth  pulse  basis  center.  Matrix  I  is  a  column  vector  which  contains  the  unknown 
coefficients  /  =  {In\  n  =  1, . . . ,  N}  that  are  used  to  approximate  the  current: 

J(p')  Pn(p')  (2.4) 

71=1 


where  Pn(p')  denotes  the  unit  pulse  basis  function  centered  at  pn.  The  column  vector  V 
elements  are  given  by  minus  the  incident  field  at  the  matching  points: 


v»  =  -E‘  (p„) 


(2.5) 


For  brevity,  the  expressions  have  been  developed  only  for  the  TM  (Ey)  polarization. 
Similar  equations  can  be  easily  derived  for  the  TE  polarization.  Nevertheless,  the  FB  method 
that  will  be  described  below  is  applicable  to  both  polarizations. 
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The  FB  method  will  be  formulated  using  the  matrix  notation  of  (2.2),  instead  of  the 
integral  equation  (2.1).  First,  consider  the  following  decomposition  applied  over  the  matrices 
involved  in  equation  (2.2): 

7  =  lf  +  f  (2.6) 

I  =  I7  + 1*  + 16  (2.7) 

where  I*  is  the  forward  component  (i.e.  the  current  contribution  due  to  the  waves  propa¬ 
gating  in  the  forward  direction) ,  Ib  is  the  backward  component  (or  current  contribution  due 

=/  =3  =b 

to  the  waves  propagating  in  the  backward  direction),  and  Z  ,  Z  and  Z  are  respectively 
the  lower  triangular  part,  the  diagonal  part  (self  impedance  terms)  and  the  upper  triangular 
part  of  Z. 

Equation  (2.2)  can  now  be  split  into  forward-propagation  and  backward-propagation 
matrix  equations,  respectively,  as  follows: 

=  F-I  f-(7f  +  Ib)  (2.8) 

T -t  =  -I6  •  (7f  +  f)  (2.9) 

where  (2.9)  has  been  assumed  by  definition.  It  can  be  seen  that,  for  a  given  nth  matching 
point  located  at  p„,  the  right  hand  side  (RHS)  of  (2.8)  contains  the  incident  field  and  the 
contribution  of  the  current  elements  located  in  the  front  of  this  receiving  element,  which 
corresponds  to  the  nth  forward  region  in  Figure  2.2(b).  Likewise,  the  RHS  of  (2.9)  contains 
the  influence  of  the  current  elements  in  the  rear  of  the  receiving  element,  so  it  represents 
the  nth  backward  region  contribution  as  in  Figure  2.2(b).  From  this  it  is  clear  that  the 
MoM  current  elements  must  be  numbered  sequentially  as  a  function  of  increasing  x  in  the 
FB  method,  i.e., 

^(Pn-l)  <  X(Pn)  <  X(Pn+ 1)'  (2‘10) 

Equations  (2.8)  and  (2.9)  can  be  solved  iteratively,  where  the  currents  in  the  ith  stage  of 
the  algorithm  (J-^,/6’^)  are  obtained  as: 


-f  ■  lm 


(2.12) 


The  algorithm  starts  with  =  0.  It  must  be  noticed  that  the  matrices  involved  in  this 
iterative  process  do  not  need  to  be  factorized  or  inverted,  because  Z  +Z*  is  a  lower  triangular 

c  h 

matrix,  and  Z  +  Z  is  an  upper  triangular  matrix;  so,  the  equations  (2.11)  and  (2.12) 
can  be  solved  for  1^  and  Ib’^  by  forward  and  backward  substitution,  respectively.  The 
convergence  has  been  shown  to  be  extremely  rapid  for  moderately  rough  surfaces,  generally 
requiring  fewer  than  10  iterations.  However,  it  should  be  noted  that  the  algorithm  may 
become  unstable  for  re-entrant  surfaces,  i.e.,  for  surfaces  where  x(pn+1)  <  x(pn)  for  one 
or  more  points  pn.  This  violates  the  sequential  number  requirement  of  (2.10).  The  GFB 
method  developed  next  overcomes  this  limitation  for  the  case  of  one  or  more  arbitrarily 
shaped  scattering  obstacles  on  the  surface. 

2.3  The  Generalized  Forward-Backward  Method 

2.3.1  Method  Description 

Consider  now  the  composite  problem  depicted  in  the  Figure  2.3,  where  one  or  more  PEC 
obstacles  (like  a  ship  or  large  rogue  breaking  wave)  are  included  in  the  surface  contour  S. 
For- this  kind  of  problem  the  conventional  FB  method  does  not  exhibit  convergent  behavior, 
because  the  presence  of  the  obstacle  highly  disturbs  the  propagation  process  assumed  by  the 
conventional  FB  method.  There  are  strong  interactions  between  the  obstacle  and  the  nearby 
ocean-like  surface,  and  within  the  obstacle  itself,  all  of  which  may  not  be  taken  into  account 
with  the  conventional  formulation  involved  in  the  standard  FB  method. 

In  order  to  overcome  this  drawback,  a  new  method  called  the  Generalized,  Forward  Back¬ 
ward  (GFB)  method  is  presented  in  this  chapter.  The  GFB  method  consists  on  a  gener¬ 
alization  of  the  conventional  FB  approach  which  enhances  the  scope  of  application  of  the 
previous  method  to  composite  scattering  problems  as  the  one  depicted  in  Figure  2.3.  This 
method  is  based  on  the  same  general  concepts  previously  stated  for  the  FB  method,  but 
includes  some  significant  differences  mainly  in  the  decomposition  of  the  matrix  Z  that  will 
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be  detailed  next.  For  the  sake  of  simplicity,  the  GFB  formulation  will  be  presented  for  a 
simple  problem  consisting  of  a  sea  surface  containing  only  one  obstacle  (see  Figure  2.3);  the 
extension  of  the  formulation  to  several  obstacles  is  obvious  as  will  be  seen  later. 

Again  starting  with  equation  (2.2),  which  has  been  obtained  from  (2.1)  after  the  MoM 
discretization  process.  In  the  same  way  as  done  in  the  FB  method,  the  current  is  expressed 
as  the  sum  of  two  contributions  ( forward  and  backward ) : 

7  =  7f  +  f  (2.13) 

but  now  the  impedance  matrix  is  split  in  a  different  way: 

f  =  l?9  +  T9  +  '  (2.14) 

where  the  Z  9  matrix  is  the  diagonal  part  of  Z  with  an  additional  block  including  the 
impedance  sub-matrix  corresponding  to  the  ship  and  near-field  sea  region  (Region  2  in  Fig¬ 
ure  2.3);  while  Z^9  and  ~Z>9  are  respectively  the  lower  triangular  part  and  the  upper  triangular 
part  of  Z  but  excluding  the  matrix  Z  9 ,  as  illustrated  in  Figure  2.3.  With  this  decompo- 
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sition,  matrix  Z  contains  both  the  self  (diagonal)  terms  and  the  interaction  of  the  whole 
obstacle  and  nearby  sea  region  together. 

Then,  the  original  system  is  transformed  in  a  similar  way  as  in  the  conventional  FB 
method,  yielding  the  following  matrix  equations: 


T-~ji  = 

(2.15) 

T-t  =  + 

(2.16) 

which  can  be  iteratively  solved  for  I^'l>  and  Ib’^  as: 

( '79+f’)l m  = 

(2.17) 

(2.18) 

starting  with  I b’^  =  0  in  (2.17). 

The  solution  of  (2.17)  and  (2.18)  differs  from  (2.11)  and  (2.12)  because  neither  Z  9  +  Z*9 

■=sg  —bg 

nor  Z  +  Z  are  triangular  matrices.  Nevertheless,  the  equations  can  also  be  easily  solved 
by  combining  forward  or  backward  substitution  together  with  the  direct  factorization  of 
the  square  block  of  Z  whose  dimension  depends  only  on  the  number  of  current  elements 
in  Region  2,  namely  the  ship  and  nearby  sea  (see  Fig.  2.3).  A  description  of  the  solution 
procedure  can  be  found  in  the  next  section.  Qualitatively,  the  solution  proceeds  as  follows. 
For  each  iteration  i: 

—f  (i) 

1.  Find  I  over  Region  1  using  the  forward  propagation  principle. 

(i) 

2.  Direct  solve  for  I  over  Region  2  using  the  fields  radiated  by  the  Region  1  currents 
plus  the  incident  field  as  excitation. 

3.  Find  I  over  Region  3  by  the  forward  propagation  principle. 

4.  Find  I6 ’  *  over  Region  3  by  the  backward  propagation  principle. 

5.  Direct  solve  for  Ib'^  over  Region  2  using  the  fields  radiated  by  the  Region  3  currents. 
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6.  Find  Ib over  Region  1  using  the  backward  propagation  principle. 

The  computational  cost  of  the  GFB  method  is  practically  the  same  as  the  conventional 
FB  method.  It  only  has  the  additional  computational  cost  of  factorizing  the  square  block 

«— ~SQ 

of  Z  which  corresponds  to  the  MoM  matrix  of  Region  2.  Nevertheless,  due  to  the  limited 
size  of  this  block  matrix,  its  factorization  can  be  performed  once  and  stored,  thus  reducing 
the  computational  work  in  subsequent  iterations  and  for  other  excitations.  So,  it  can  be 
concluded  that  the  GFB  method  has  a  computational  cost  of  0(N2)  per  iteration  as  in 
the  FB  method  [6].  The  storage  requirement  is  0(N )  to  store  the  iterated  currents  as  in 
the  FB  method,  with  the  additional  storage  of  a  square  matrix  of  size  M  x  M,  where  M 
is  the  number  of  current  elements  included  in  the  MoM  region  (Region  2  in  Figure  2.3). 
Generally,  M  <C  N,  but  the  0(M 2)  matrix  storage  requirement  may  be  comparable  to  the 
O(N)  storage  of  the  currents.  It  is  noted  that  the  0(N )  storage  requirement  assumes  that 
the  full  Z  matrix  is  not  stored,  so  it  is  necessary  to  recompute  the  matrix  elements  at  each 
iteration  (except  for  the  inner  MoM  matrix  associated  with  Region  2).  However,  the  simple 
closed  form  of  the  matrix  elements  allows  them  to  be  computed  very  quickly,  and  the  number 
of  iterations  is  usually  very  small  (typically  less  than  10).  Finally,  it  is  noted  that  the  new 
spectral  technique  introduced  in  [6]  for  accelerating  the  FB  method  may  likewise  be  used  to 
accelerate  the  GFB  method,  with  some  minor  modifications  as  described  in  Section  2.4. 

2.3.2  Solving  equations  (2.17)  and  (2.18) 

In  this  section  a  description  of  the  procedure  used  to  solve  equations  (2.17)  and  (2.18)  is 

~sg  =fg  =sg  =bg 

presented.  It  must  be  pointed  that  neither  Z  +  Z  nor  Z  +  Z  are  triangular  matrices; 
nevertheless,  these  equations  can  be  easily  solved  by  combining  forward  or  backward  substi¬ 
tution  together  with  the  direct  factorization  of  a  square  block  sub-matrix,  whose  dimension 
is  just  the  extent  of  the  MoM  region.  For  brevity,  we  will  explain  only  the  resolution  of 

=sg  ~fg 

(2.17)  which  concerns  the  matrix  Z  +Z  .  Equation  (2.18)  can  be  solved  in  a  similar  way. 
Equation  (2.17)  has  the  general  form 

2-x  =  b  (2.19) 
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Figure  2.4:  Block  decomposition  of  a  “quasi”  lower  triangular  matrix. 


where  A  =  z'9  +  Zf9 ,  x  =  7/,W,  and  b  =  V  -  Z?9  •  A  is  a  “quasi”  lower  triangular 

matrix  as  depicted  in  Figure  2.4.  In  order  to  solve  these  equations  in  an  efficient  way,  the 
matrices  A  ,x  and  b  of  (2.19)  are  subdivided  into  blocks  as  shown  in  Figure  2.4.  Thus  (2.19) 
can  be  expressed  in  terms  of  the  blocks  as  follows: 

(%i+Xi)'Xi  =  bi  (2.20) 

A21  '  X\  +  A22  •  (2.21) 

^4-31  ■  +  A32  •  X2  +  ^^33  +  ^33^  •  %3  =  ^3  (2.22) 

The  solution  is  then  obtained  as  follows: 

1.  Equation  (2.20)  is  a  lower  triangular  matrix,  so  X\  can  be  easily  obtained  by  forward- 
substitution. 


2.  Once  Xi  is  has  been  calculated,  X2  can  be  obtained  from  (2.21)  by  the  direct  solution 


A22  •  x2  ~  {i>2  ~  A21  ‘  Xi ) 


(2.23) 


LU  decomposition  may  be  used  so  that  A22  needs  to  be  factorized  only  once  for  a  given 
Region  2  geometry.  The  factorized  matrix  may  then  be  saved  for  subsequent  iterations 
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and  for  other  excitations. 


3.  Finally,  once  x\  and  X2  have  been  calculated,  x$  can  be  obtained  from  (2.22)  again  by 
forward  substitution  from: 

f-A33  +  A33^  ’  —b^  —  A$\  •  £1  —  A32  •  %2  (2.24) 

- f 

because  A33  +  A33  is  a  lower  triangular  matrix. 

A  similar  procedure  is  applied  for  the  solution  of  equation  (2.18).  The  only  difference  is 
that  in  this  case  we  will  define  upper  triangular  component  matrices  that  can  be  solved  by 
back  substitution. 

2.3.3  Numerical  results 

In  this  section  results  are  presented  to  validate  the  convergence  and  accuracy  of  the  new  GFB 
method,  and  to  investigate  some  of  the  effects  of  a  random  rough  surface  on  the  backscatter 
pattern  of  a  target  on  the  surface.  In  the  following,  the  sea  surfaces  are  randomly  generated 
realizations  of  a  Gaussian  random  process  with  a  Pierson-Moskowitz  spectrum  for  a  given 
wind  speed  [8].  The  excitation  antenna  is  a  linear  array  of  15  equally  spaced  electric  line 
sources  with  a  cosine  aperture  weighting  function.  This  source  array  produces  a  9°  main 
beam  with  low  sidelobes.  (Low  sidelobes  are  important  so  that  the  sea  surface  near  the 
source  is  not  strongly  illuminated.)  The  polarization  is  horizontal,  as  in  the  derivation  for 
the  TM  EFIE  earlier. 

Figure  2.5  shows  the  source  location  and  pattern  function  relative  to  a  PEC  ship-like 
target  on  a  PEC  rough  sea  surface.  For  this  example,  the  sea  surface  is  204.8  in  in  length, 
the  wind  speed  is  15  meters  per  second  (m/s),  and  the  wavelength  is  1  m.  The  source  is  10  m 
high  on  the  extreme  left,  and  points  horizontally.  The  MoM  region  shown  within  the  dashed 
vertical  lines  contains  the  ship  and  5  m  of  sea  surface  on  either  side.  For  validating  the  new 
method  with  this  example,  the  current  obtained  with  GFB  is  compared  with  the  solution 
given  by  MoM  applied  to  the  whole  surface  in  Figure  2.6,  as  a  function  of  the  pulse  basis 
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Figure  2.5:  Source  and  target  geometry  on  a  sea  surface  with  wind  speed  15  m/s.  Units  are 
in  meters. 


index  number.  For  the  GFB  results,  the  MoM  region  is  within  the  outer  dashed  vertical 
lines,  and  the  ship  is  within  the  inner  dashed  lines.  10  iterations  are  used.  The  difference 
between  the  two  currents  is  shown  in  the  Figure  2.7.  After  only  10  iterations,  the  maximum 
difference  in  the  currents  is  on  the  order  of  10-5. 

The  residual  error  is  used  for  monitoring  the  convergence  of  the  GFB  in  terms  of  the 
number  of  iterations.  The  residual  error  vector  after  the  ith  iteration  is  defined  as: 

rw  =  F-I-7W.  (2.25) 


By  substituting  equations  (17)  and  (18)  in  (2.25),  the  residual  error  vector  can  be  evaluated 
in  a  more  efficient  way  as: 


f(i)  _  g19  .  [yh0-i)  _  jKO 


(2.26) 


The  Residual  Error  is  defined  as: 

f(0 

Residual  Error  =  _ 

V 


(2.27) 


where  |  r^  denotes  the  vector  norm.  The  residual  error  of  the  GFB  for  this  example  is 
shown  in  Figure  2.8;  it  decreases  exponentially  with  the  number  of  iterations.  The  residual 
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Current  magnitude  Current  magnitude 


Figure  2.6:  Magnitude  of  the  currents  as  a  function  of  index  number  (with  10  segments/A): 
(a)  the  GFB  Method  (10i/l  iteration)  and  (b)  Moment  Method. 
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Figure  2.7:  Magnitude  of  the  difference  between  GFB  and  MoM  currents  as  a  function 
index  number. 
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Figure  2.8:  Residual  and  absolute  errors  versus  the  number  of  iterations. 


error  of  the  FB  method  has  also  been  plotted,  and  as  for  reasons  explained  in  Section  3,  the 
FB  method  generally  cannot  achieve  a  convergent  solution  when  an  obstacle  is  present  in 
the  sea  surface.  Figure  2.8  also  plots  the  Absolute  Error  of  the  GFB  method,  defined  by 


Absolute  Error  = 


7(i 

—  I  MoM 

I  MoM 

(2.28) 


where  ImoM  is  the  MoM  reference  solution  for  the  currents.  It  is  seen  that  the  absolute 
error  also  falls  exponentially,  but  after  about  14  iterations  it  levels  out.  This  is  often  a 
characteristic  of  iterative  methods  in  general,  and  is  caused  by  the  finite  numerical  precision 
of  the  computer.  To  judge  convergence  of  the  solution,  the  residual  error  is  a  very  good 
indicator  and  is  easy  to  compute  at  each  iteration  using  (2.26)  and  (2.27).  Halting  the 
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algorithm  when  the  residual  error  reaches  about  10~2  or  10-3  yields  very  accurate  scattering 
results.  For  this  example,  6-10  iterations  is  found  to  be  quite  sufficient. 

To  investigate  the  scattering  behavior  of  a  ship  on  a  rough  surface,  the  geometry  of 
Figure  2.9  is  considered.  The  ship  and  the  source  are  the  same  as  in  Figure  2.5,  but  the 
length  of  the  sea  surface  is  409.6  m,  and  the  source  antenna  moves  in  elevation  at  a  constant 
distance  from  the  ship  to  generate  a  backscatter  pattern  (i.e. ,  the  pattern  of  the  scattered 
field  at  the  source  position  as  a  function  of  elevation  angle).  The  antenna  beam  points 
directly  at  the  ship,  and  the  source  starts  at  a  height  of  10  m  above  the  surface.  The 
wavelength  is  again  1  m.  Two  different  randomly  generated  sea  surfaces  at  the  same  wind 
speed  (15  m/s)  are  shown  in  the  figure.  It  is  noted  that  the  ship  remains  upright  and  does 
not  roll  with  the  waves  in  the  results  that  follow. 


Figure  2.9:  Geometry  for  backscatter  patterns.  In  the  pattern  plots,  continuous  line  repre¬ 
sents  Surface  #1  and  dotted  line  Surface  #2.  (The  ship  is  the  same  in  both  cases.)  Wind 
speed  =  15  m/s. 
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Figure  2.10:  Comparison  between  the  backscattered  field  of  a  ship  on  a  rough  surface  (Sur¬ 
face  #1)  and  on  a  flat  surface.  A  reference  solution  for  an  infinite  flat  surface  found  using 
MoM  and  image  theory  is  also  included. 

Figure  2.10  shows  the  backscatter  patterns  for  the  ship  on  Surface  #1,  and  for  the  ship 
on  a  finite  flat  surface  of  the  same  length  for  comparison.  Also  shown  is  a  reference  MoM 
solution  for  the  ship  on  an  infinite  flat  surface  found  using  image  theory  [4].  The  very  close 
agreement  between  the  image  theory  result  and  the  GFB  result  for  the  finite  flat  surface  show 
that  end  point  effects  are  negligible,  and  further  validates  the  GFB  method  for  scattering 
problems.  The  figure  also  shows  that  the  rough  surface  greatly  affects  the  backscatter  pattern 
compared  with  a  flat  surface. 
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0  (deg) 

Figure  2.11:  Comparison  between  the  backscattered  field  for  Surfaces  #1  and  #2  with  the 
ship.  The  backscattering  from  surface  #1  without  the  ship  is  also  plotted. 

Figure  2.11  shows  the  backscatter  patterns  for  the  ship  on  the  two  rough  surfaces  depicted 
in  Figure  2.9.  It  is  seen  that  the  patterns  are  significantly  different  for  these  two  surfaces 
which  have  the  same  roughness  scale,  although  it  could  be  argued  that  the  patterns  are  at 
about  the  same  average  level.  The  backscatter  due  to  Surface  #1  without  the  ship  is  also 
plotted  to  show  the  level  of  the  background  clutter  for  this  wind  speed.  The  clutter  level 
computed  here  is  quite  low,  but  it  should  be  mentioned  that  this  ocean  model  does  not 
include  the  effects  of  breaking  waves  which  generally  dominate  the  backscatter  return  from 
the  ocean  in  the  absence  of  a  ship.  In  fact,  the  clutter  level  in  Figure  2.11  is  due  almost 
entirely  to  direct  reflection  of  the  downward  pointing  low  sidelobe  beams.  However,  as  noted 
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Figure  2.12:  Comparison  between  the  backscattered  field  for  a  ship  on  Surface  #1  and  on  a 
randomly  generated  surface  with  wind  speed  =  7  m/s. 

earlier,  the  GFB  method  could  be  used  to  analyze  an  isolated  breaking  wave  in  the  presence 
of  an  otherwise  smoothly  varying  ocean  surface  simply  by  replacing  the  ship  geometry  with 
that  of  the  wave. 

Figure  2.12  shows  the  backscatter  patterns  for  the  ship  on  Surface  #1  and  on  a  randomly 
generated  surface  with  wind  speed  7  m/s.  Again,  the  patterns  are  significantly  different, 
although  they  are  at  about  the  same  average  level.  Furthermore,  the  7  m/s  result  does  not 
appear  any  closer  to  the  flat  surface  result  of  Figure  2.10  than  the  15  m/s  result,  even  though 
the  roughness  scale  is  smaller.  Any  further  conclusions  about  these  results  requires  a  more 
thorough  statistical  study,  which  is  the  subject  of  Chapter  5. : 
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2.3.4  Discussion  and  Conclusions 

The  new  approach  has  been  shown  to  provide  very  accurate  results  and  maintains  the  same 
fast  convergence  (usually  in  less  than  10  iterations)  and  0(N2)  computational  cost  associ¬ 
ated  with  the  conventional  FB  solution.  The  storage  requirement  is  0(N  +  M2),  where  M 
is  the  number  of  basis  functions  used  in  the  small  MoM  region  around  the  target.  Neverthe¬ 
less,  the  0(N 2)  computational  cost  makes  the  method  inefficient  for  very  large  surfaces.  In 
order  to  overcome  this  limitation  the  GFB  could  be  combined  with  modern  integral  equa¬ 
tion  acceleration  algorithms  such  as  the  asymptotic  fast  multipole  method  [9],  or  the  novel 
spectral  acceleration  method  [6]  which  could  reduce  the  computational  cost  to  O(N).  The 
improvement  of  the  GFB  method,  combining  it  with  novel  acceleration  algorithms,  will  be 
addressed  in  the  next  section. 

The  numerical  results  show  that  different  random  ocean  surfaces  greatly  affect  the  backscat- 
ter  patterns  of  a  ship-like  target,  although  the  patterns  tend  to  have  about  the  same  average 
level.  It  is  difficult  to  make  any  further  conclusions  at  this  point  about  how  the  wind  speed 
affects  the  backscatter  patterns.  One  would  expect,  for  example,  that  higher  wind  speeds 
(and  therefore,  larger  roughness  scales)  would  cause  a  larger  variance  in  the  backscatter.  A 
more  thorough  statistical  investigation  is  necessary,  which  is  the  subject  of  Chapter  5.  The 
purpose  of  the  GFB  method  introduced  here  is  to  provide  a  numerical  tool  for  studying 
this  class  of  scattering  problems,  and  to  provide  a  reference  solution  for  more  approximate 
techniques  discussed  later  in  this  report. 

2.4  Spectral  Acceleration  for  the  GFB 

In  the  Generalized  Forward-Backward  method,  the  direct  computation  of  the  matrix-vector 
products,  included  in  (2.17)  and  (2.18),  Z*9  •  (J^  +  /6,(l  ^ )  and  Z  9  •  (1^  +1  ’^)  implies 
a  computational  cost  of  0(N2).  In  [6]  a  spectral  acceleration  for  the  FB  was  introduced. 
This  acceleration  reduces  the  computational  cost  of  these  products  to  O(N),  and  is  modified 
here  in  order  to  apply  it  in  the  GFB  method.  In  this  section  the  modifications  are  described 
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and,  for  clarity,  the  description  is  divided  into  3  parts  which  correspond  to  the  3  regions  in 
Figure  2.3.  For  brevity,  only  the  acceleration  for  the  forward  propagation  is  included;  the 
formulation  for  the  backward  propagation  is  obtained  in  the  same  way. 

2.4.1  Implementation  of  the  Spectral  Acceleration 

Region  #1 

In  order  to  evaluate  1^  in  (2.17)  at  the  points  in  Region  #1  (n  <  Nx),  a  matrix-vector 
product  of  the  form 

E/W.)  =  E  where  Z^  =  H^(k\p„  -  pm\)  (2.29) 

m=l 

must  be  computed.  In  the  above  expression,  Im  =  4- (It  is  noted  that  the 

constants  relating  Hq  to  the  2D  Green’s  function  are  absorbed  into  Im). 

To  reduce  the  computational  cost  of  this  computation,  the  spectral  acceleration  presented 
in  [6]  divides  the  elements  into  two  groups  for  a  given  receiving  element:  strong  interaction 
and  weak  interaction  groups,  as  shown  in  Figure  2.13(a).  The  criterion  that  defines  these 
groups  is  the  distance  from  the  receiving  element.  For  a  given  distance  Ls  which  includes 
Ns  elements  (Ls  =  Ns  •  Ax ,  where  Ax  is  the  sampling  increment),  the  forward  contribution 
is  decomposed  as: 

E^(pn)  =  Ew(pn)  +  Es(pn) 

E*(Pn)  = "  £  ‘  In  Zll  ;  E,(M  =  £  Im  Z&  (2.30) 

m= 1  m~n—Ns 

where  Ew  and  Es  are  the  fields  of  the  weak  and  strong  interaction  groups,  respectively. 

The  radiation  of  the  strong  interaction  group  is  computed  as  in  (2.30),  but  the  weak  group 
contribution  is  obtained  by  employing  the  spectral  representation  of  the  Hankel  function: 

Hi2)(k\pn  -  pm\ )  =  ~  [  e-Mxn-xm) cos 4>+(zn-zm) sin*]  d(j)  =  (2.31) 

7T  JC,j, 

where  C <p  is  the  path  of  integration  in  the  complex  <f)  plane  (shown  later  in  Figure  2.14). 
Introducing  (2.31)  into  the  expression  for  Ew(pn )  and  interchanging  the  integration  and 
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summation,  yields, 


(2.32) 


K(M  =  -[  FM 

7T  J  C# 


where 


e~jk[(xn-xm )  cos  4>-Zm  sin  <j>\ 


(2.33) 


It  is  at  this  point  where  the  reduction  in  the  computational  cost  is  obtained  because  it  can 
be  shown  that  Fn(</>)  is  easily  updated  using  a  recursive  procedure: 


FJcp)  =  Fn-M  e-M*n-xn-i)cos<t, 


Up  to  this  point  the  spectral  acceleration  procedure  is  the  same  as  presented  in  [6].  The 
modifications  will  appear  in  regions  #2  and  #3  as  discussed  next,  and  in  the  numerical 
evaluation  of  the  integral  in  (2.32)  which  will  be  discussed  at  the  end  of  this  section. 

Region  #2 

Referring  to  Figure  2.13.b,  the  sampled  values  for  the  current  in  this  region  (J^’W,  Ni  <  n  < 
N\  +  N2)  are  obtained  with  the  direct  solving  of  the  following  matrix  system: 

Xl  ■  7n(i)  =  V~  Ef  n,me  (JVx , N,  +  N2]  (2.35) 

The  evaluation  the  right-side  of  (2.35)  requires  the  computation  of  E^\pn)  for  each  of 
the  elements  of  this  region  using, 

JVi 

E?(Pn)  =  ZS..  (2.36) 

m=  1 

Due  to  the  different  way  in  which  the  current  is  evaluated  in  this  region  (compared  with 
regions  #1  and  #3),  the  group  division  is  done  using  another  criterion.  An  element  belongs 
to  the  weak  group  if  the  distance  from  this  element  to  the  first  element  in  region  #2  (element 
#(IVi  +  1))  is  larger  than  Ls.  So,  in  terms  of  number  of  elements,  it  can  be  express  as: 

Ni-N,  Nx 

E(;\M=  Y,  E  1™Z!X  =  EV(M  +  E?(i>j  (2.37) 

m=  1  m=N1-Ar5  +  l 
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It  must  be  noticed  that  with  this  criterion,  the  weak  and  the  strong  groups  remain  constant 
for  all  the  elements  of  the  region  #2.  Taking  this  point  into  account,  the  recursive  procedure 
to  update  the  function  Fn(4>)  must  be  changed  to: 

Fn(4>)  =  Fn-i(<f>)  n  £  (Ni,Ni  +  jV2]  (2.38) 

Region  #3 

In  the  region  #3  (n  >  N\  +  N2),  the  forward  current  component  is  obtained  by  forward 
substitution  as  in  the  region  #1  and  the  following  matrix-vector  product  must  be  computed: 

E,(M  =  E  Zll-  (2-39) 

771“  1 

Therefore  the  procedure  used  to  evaluated  the  contribution  of  the  weak  group  elements  is 
going  to  be  similar  to  that  in  Region  #1.  Some  modifications  are  included  for  the  first  Ns 
elements  of  this  region  for  updating  the  weak  group  due  to  the  presence  of  the  MoM/ship 
region  (Region  #2).  This  splits  the  Region  #3  computation  into  two  cases  shown  in  Figures 
2.13  (c)  and  (d). 

(1)  Figure  2.13(c):  n  <  Ni  +  N2  +  Ns 

'  For  the  first  Ns  elements  of  this  region  a  special  division  for  the  strong  group  is  con¬ 
sidered.  The  strong  group  will  be  formed  by  all  the  elements  of  Region  #2,  plus  the 
Ns  elements  in  Regions  #1  and  #3  which  precede  the  receiving  element.  Therefore, 
the  strong  group  will  contain  N2  +  Ns  elements.  With  this  distribution,  (2.39)  can  be 
rewritten  as: 

7i-l— N2—Ns  n—  1 

E,(Pn)=  E  ImZl’n+  E  Zfa  =  EW(M  +  E,(Pn)  (2.40) 

771=1  m~n—N2-Ns 

and  the  procedure  to  obtain  F„(</>)  must  change  to: 

Fn((f))  =  Fn-l(<l>)  e-Mxn-Xn-l)COS<p  +  .  e~jk(xn-X  ^  )  COS  <l>  ^jkzj  sin  *  (2.41) 

where  n  =  n  —  1  —  Ns  —  N2- 
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(2)  Figure  2.13(d):  n  >  N\  +  N2  +  Ns 


When  a  receiving  point  if  this  set  is  considered,  all  the  elements  of  the  strong  group  are 
placed  in  region  #3,  and  the  elements  of  region  #2  will  be  a  part  of  the  weak  group. 
Therefore,  (2.39)  can  be  divided  as  in  (2.30)  as 

=  (2.42) 

m= 1  m=n-Ns 

To  compute  Ew(pn),  the  function  Fn{(j))  can  be  evaluated  as  in  (2.34),  after  updating 
its  value  at  the  first  element  of  this  set  (p  =  Ni  +  N2  +  Ns  + 1)  with  the  contributions 
due  to  the  elements  in  region  #2: 

N!+N2 

Fp(4>)  =  Fp_i(0)  4'  •  e-jkiXp-x»)cos4,  e"J'feVsin*  (2.43) 

n  1 

For  the  following  points  (n  >  p),  the  values  of  Fn((f))  are  obtained  via  (2.34). 

2.4.2  Contour  Integral 

The  previous  points  of  this  section  were  focused  on  how  to  evaluate  the  integrand  of  (2.32) 
for  the  different  regions.  In  this  part,  the  numerical  computation  of  this  integral  is  explained, 
including  how  the  integration  path,  the  limits  and  the  numerical  sampling  density  are  de¬ 
termined.  The  procedure  to  obtain  these  three  parameters  is  basically  the  same  used  in  [6], 
although  some  modifications  are  included  due  to  the  presence  of  a  target  on  the  sea  surface. 
These  modifications  affect  mainly  the  integration  path  and  the  limits  of  integration. 

As  explained  in  [6],  because  of  the  Hankel  function  is  analytic  in  the  complex  space, 
the  integration  contour  C $  can  be  deformed  to  another  integration  path  C$  as  shown  in 
Figure  2.14.  This  path  is  chosen  in  order  to  reduce  the  computational  cost  needed  to  evaluate 
the  integral,  and  to  reduce  errors  caused  by  the  limited  numerical  precision  of  the  computer 
when  the  integrand  presents  large  values.  Cs  follows  a  straight  line  in  the  range  —<f)s,max  < 
Re{(j)}  <  (j>s,max  where  <4, max  is  defined  by 

,  ,  -l  /  zmax  zmin\  fo  A  A\ 

4, max  =  tan  - - - J  ,  (2.44) 
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and  is  illustrated  in  Figure  2.15.  (Notice  that  <fis,max  is  the  maximum  angle  from  horizontal 
formed  by  any  point  of  the  weak  group  and  the  receiving  element).  S  is  the  angle  that  Cs 
makes  with  the  real  axis  of  the  complex  (f>  plane,  and  is  chosen  such  that  the  integrand  of 
(2.32)  does  not  become  too  large  anywhere  along  Cs- 

The  steepest  descent  path  (SDP)  through  the  saddle  point  4>s,max  is  also  in  the  Figure  2.14, 
along  with  the  steepest  ascent  path  which  is  orthogonal  to  the  SDP  at  the  saddle  point.  It 
is  noted  that  the  SDP  and  SAP  make  an  angle  of  45°  with  the  real  axis.  For  different  sets 
of  source  and  receiver  elements,  the  saddle  point  moves  to  different  points  along  the  Re{<f>} 
axis  in  the  range  —<j)s,max  <  Re{(j>}  <  <t>s,max-  The  integrand  will  have  its  maximum  value  on 
Cs  at  the  point  where  the  SAP  through  (f>s,max  crosses  Cs  (point  A  in  Figure  2. 14). 

The  angle  5  is  given  by  AOB,  where  B  =  <j)s,max  and  A  is  the  point  on  the  SAP  of  (j)s,max 
where  the  integrand  is  equal  to  e2.  This  criterion  guarantees  that  the  every  contribution 
of  the  elements  in  the  weak  group  are  smaller  than  e2,  which  is  chosen  to  keep  the  integral 
numerically  stable.  It  is  noted  the  in  [6]  this  limiting  value  was  chosen  to  be  e20.  The  reason 
a  much  lower  is  used  here  is  because  many  of  the  saddle  points  associated  with  the  Region 
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#2  target  structure  are  grouped  close  to  </>S)mox  (see  Figure  2.15),  so  their  total  contribution 
to  the  integrand  via  (2.33)  will  be  much  larger  than  e2.  In  [6]  the  saddle  points  are  grouped 
more  around  the  origin,  with  only  a  few  point  near  ±(f)s,max ■  With  this  new  limit  of  e2,  it 
may  be  shown  that  8  is  given  approximately  by, 


5  —  tan  ^l/fc) 


with 


1 


(2.45) 


where  Rs  =  Jl%  +  (zmax  +  zmin)2.  Using  this  new  limit,  the  integrand  in  (2.32)  does  not 

MoM 


present  very  large  values  as  demonstrated  in  Figure  2.16  for  a  receiving  position  at  the  top 
of  the  target.  This  allows  the  integral  to  be  evaluated  very  accurately  as  is  shown  in  the 
numerical  results. 

Once  the  integration  path  has  been  determined,  the  next  step  is  to  set  up  the  integration 
limits.  These  limits  can  be  easily  obtained  by  using  geometrical  considerations.  If  </>s>max 
is  defined  as  in  Figure  2.15,  where  zmax  and  zmin  will  denote  the  maximum  and  minimum 
z  coordinates  of  the  whole  surface  including  the  target,  then  the  integration  limits  will  be 
[— <f>s,max,  <f>s,max\-  Figure  2.16  plots  the  integrand  of  (2.32)  for  a  point  at  the  top  of  the 
target,  it  can  be  noticed  that  the  main  contributions  are  inside  the  integration  interval,  and 
the  integrand  is  nearly  zero  outside. 

The  value  of  the  integration  step  size  used  for  the  FB  method  in  [6]  is  also  used  here, 
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Figure  2.16:  Integrand  of  (10)  for  a  point  at  the  top  of  the  target. 


given  by 

A*  =  g/22  (2.46) 

and  the  variable  of  the  integration  is  mapped  to  the  Recj)  axis  according  to 

d<f>  A(j)ejS  (2.47) 


Using  the  above  integration  limits  and  step  size,  the  integral  of  (2.32)  it  is  computed  by 
adding  the  contribution  of  2  Q  + 1  plane  waves  in  which  the  integral  is  discretized.  The  value 
of  Q  is  obtained  from: 


Q  = 


A 

Ac j> 


where  [•]  denotes  the  nearest  integer  larger  than  the  argument. 


(2.48) 


2.4.3  Computational  cost 

The  computational  cost  associated  to  the  Generalized  Forward- Backward  Method  is  0(N 2) 
per  iteration.  This  includes  the  0(N 2)  forward  (or  backward)  integration  of  the  currents 
radiated  to  every  element,  plus  the  0(1V|)  computational  cost  of  a  matrix- vector  product  to 
obtain  the  current  over  the  MoM  Region  #2.  (The  Region  #2  MoM  matrix  is  assumed  to 
be  factorized  prior  to  the  iteration  process.)  Considering  that  N2  N,  the  computational 
cost  per  iteration  of  the  Generalized  Forward-Backward  method  is  0(N 2) 
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When  the  Spectral  Acceleration  method  described  in  this  section  is  applied  to  the  GFB, 
the  computational  cost  per  forward  or  backward  iteration  is  estimated  as  follows.  The  cost 
for  the  strong  group  iteration  for  N  receiving  points  is  Ns  •  N.  (One  unit  of  cost  is  defined 
roughly  as  one  complex  exponentiation  plus  one  complex  multiplication  plus  one  complex 
addition.)  The  cost  of  computing  the  spectral  amplitudes  (2.34)  of  2 Q  +  1  plane  waves  is 
2(2Q  + 1)  •  (N  —  Ns),  and  the  cost  of  computing  the  integral  (2.32)  is  (2 Q  + 1)  •  (N  —  Ns)  for 
N  -  Ns  receiving  elements.  The  cost  of  the  Region  #2  matrix  vector  product  is  again  N$. 
The  total  cost  per  iteration  is  then  given  approximately  by 

Cost  «  N  •  Ns  + 3  (2Q  +  1)  •  (iV  —  Ns)  4-  JVf .  (2.49) 

Since  7V2,  Ns,  and  Q  are  fixed  for  a  given  target  size,  frequency,  and  surface  roughness,  the 
total  cost  scales  as  O(N)  as  the  length  of  the  surface  increases  (provided  N$  N  '  Ns). 

As  a  result  of  the  introduction  of  the  spectral  acceleration,  the  computational  cost  of  the 
GFB  is  reduced  from  0(N2)  to  O(N)  when  increasing  surface  lengths  for  a  given  frequency 
are  considered.  For  the  case  of  a  given  target  and  a  given  sea-surface,  the  computational 
cost  as  a  function  of  frequency  is  0(N 2)  because  Ns  scales  with  N  for  this  case.  However, 
even  in  this  situation  the  spectral  acceleration  greatly  reduces  the  computational  cost  of  the 
method.  Compared  with  the  ordinary  GFB,  for  a  given  problem  size  the  accelerated  version 
(fast  GFB)  reduces  the  computational  cost  by  a  factor  of  about  Ls/L  where  L  is  the  total 
length  of  the  surface 

2.4.4  Results 

Results  are  included  in  this  section  to  demonstrate  the  accuracy  and  efficiency  of  the  spectral 
acceleration  introduced  into  the  GFB  method  for  different  incident  fields,  frequencies  and 
target  heights.  For  all  the  results,  a  Pierson-Moskowitz  ocean  spectrum  wave  model  is  used 
to  generate  random  sea  surfaces. 

Figure  2.17  shows  the  magnitude  of  the  computed  currents  as  a  function  of  the  index 
number  for  a  13  m  high  target  on  a  409.6  m  sea  surface  with  h  wind  speed  of  15  m/s  (RMS 
wave  height  =  1.215  m).  The  incident  field  is  produced  by  a  linear  array  of  15  equally  spaced 
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electric  line  sources  with  a  cosine  aperture  weighting  function  to  reduce  the  sidelobe  levels. 
The  source,  placed  at  344  m  from  the  center  of  the  ship,  points  directly  to  the  ship  and 
forms  a  5°  angle  with  the  horizontal  axis.  The  geometry  and  source  location  are  drawn  to 
scale  in  Figure  2.17(a).  Figure  2.17(b)  shows  the  currents  at  300  MHz  computed  using  the 
fast  GFB  method,  and  Figure  2.17(c)  shows  the  difference  between  these  currents  and  the 
currents  computed  by  the  original  GFB  method.  Figures  2.17(d)  and  (e)  show  the  same 
results  as  (b)  and  (c),  respectively,  at  1.2  GHz.  Figure  2.18  shows  the  same  results  at  300 
MHz  for  a  target  that  is  a  26  m  high. 

It  is  clear  that  the  relative  error  introduced  by  the  spectral  acceleration  for  these  cases 
is  on  the  order  of  10-3  or  less.  This  shows  that  the  combination  of  the  GFB  method  with 
spectral  acceleration  provides  an  accurate  method  to  analyze  the  scattering  from  targets  on 
very  large  sea  surfaces. 

In  the  previous  examples,  the  incident  field  was  generated  by  an  array  whose  radiation 
pattern  was  designed  to  minimize  the  diffraction  effects  at  the  endpoints  of  the  finitely 
truncated  sea  surface.  In  the  following  examples,  an  exponentially  tapered  resistive  sheet 
(R-card)  will  be  placed  at  both  ends  of  the  sea  surface  to  minimize  the  end  point  effects  as 
done  in  [10].  This  allows  a  pure  plane  wave  incident  field  to  be  used  instead  of  a  localized 
directive  source. 

With  the  R-cards  present,  the  surface  impedance  near  the  ends  varies  as 

A ZmTn  =  Ri  ■  (eQX  -  1)  where  a  =  (2.50) 

x  is  the  distance  from  the  junction  between  the  R-card  and  the  sea  surface,  and  Lr  is  the 
length  of  the  R-card.  The  surface  impedance  is  zero  where  the  R-card  joints  the  PEC  sea 
surface,  and  increases  exponentially  to  R2  —  R\  at  the  ends.  The  values  of  Ri  and  R2  used 
in  the  results  which  follow  are 

Ri  =  58.5  x  10"°-49  Q/ □ 

R2  =  58.5  x  10°-55  n/U  (2.51) 

The  surface  impedance  is  easily  incorporated  into  the  impedance  matrix  simply  by  adding 


36 


0  500  1000  1500  2000  2500  3000  3500  4000 

Index  number 


0  2000  4000  6000  8000  10000  12000  14000  16000 

Index  number 


Index  number 

Figure  2.17:  (a)  Example  geometry  and  source  position  (units  in  meters)  for  a  13  m  high  target,  (b) 
Magnitude  of  the  current  in  A/m  obtained  by  applying  the  Spectral  Acceleration  (10t/l  iteration)  at 
300  MHz,  and  (c)  Magnitude  of  the  difference  between  the  previous  result  and  a  reference  solution 
(GFB)  (10t/l  iteration),  (d)  Same  result  as  (b)  at  1.2  GHz.  (e)  Same  result  as  (c)  at  1.2  GHz.  The 
length  for  the  strong  group  is  Ls  =  12 m. 
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Figure  2.18:  (a)  Example  geometry  and  source  position  (units  in  meters)  for  a  26  m  high  target,  (b) 
Magnitude  of  the  current  in  A/m  obtained  by  applying  the  Spectral  Acceleration  ( 1 0th  iteration)  at 
300  MHz,  and  (c)  Magnitude  of  the  difference  between  the  previous  result  and  a  reference  solution 
(GFB)  (10th  iteration).  The  length  for  the  strong  group  is  Ls  =  24  m. 
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Figure  2.19:  Distribution  of  the  currents  with  and  without  the  r-cards. 


Zs  to  the  self  terms  Zmm  in  the  R-card  regions. 

Figure  2.19  shows  the  currents  computed  with  and  without  the  r-cards  for  a  plane  wave 
illuminating  the  target  of  Figure  2.17(a)  on  a  409.6  m  flat  surface.  The  elevation  angle  of 
the  plane  wave  is  5°  and  its  frequency  is  300  MHz.  The  R-card  placed  the  left  end  of  the  sea 
surface  is  50 A  long  and  the  one  placed  at  the  right  end  is  10A.  It  is  seen  that  the  R-cards 
minimize  the  end-point  effect  by  smoothly  tapering  the  current  to  zero.  The  effect  of  the 
R-cards  on  the  scattered  fields  can  be  seen  in  Figure  2.20  where  the  bistatic  2D  radar  cross 
section  (RCS)  due  to  both  currents  is  plotted.  The  use  of  the  R-cards  removes  the  rapid 
variations  caused  by  the  endpoint  diffractions. 

One  of  the  reasons  to  include  the  spectral  acceleration  in  the  GFB  method  is  to  be 
able  to  analyze  very  large  sea  surfaces  when  the  incidence  angle  is  very  low,  because  grazing 
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Figure  2.20:  Bistatic  RCS  due  to  the  currents  with  and  without  the  r-cards.  Frequency  =  300  MHz, 
elevation  angle  =  5°. 

incidence  angles  tend  to  illuminate  a  very  large  portion  of  the  sea  surface.  The  next  examples 
show  how  the  truncated  length  of  the  sea  surface  affects  the  RCS  computations. 

Figure  2.21  shows  the  bistatic  RCS  for  the  same  flat  surface  case  as  Figure  2.20  for  two 
finite  surface  lengths:  409.6  m,  3276.8  m.  The  lengths  of  the  surfaces  in  front  of  the  target 
are  340m  and  3200  m,  respectively.  Again,  the  plane  wave  elevation  angle  is  5°  and  the 
frequency  is  300  MHz.  The  infinite  surface  result  is  found  using  image  theory.  Both  finite 
surface  results  match  nearly  exactly  the  image  theory  reference  solution,  but  it  is  clear  that 
the  larger  surface  is  more  accurate  at  very  low  angles.  The  numerical  precision  of  the  fast 
method  and  the  tapered  R-card  endpoint  treatment  is  quite  impressive  considering  that  the 
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Figure  2.21:  Bistatic  RCS  for  two  different  lengths  of  a  flat  surface.  Frequency  =  300  MHz, 
elevation  angle  =  5°. 

currents  over  target  plus  the  entire  sea  surface  are  integrated  to  give  the  scattered  field. 

In  the  Figure  2.22,  the  same  problem  is  analyzed  but  in  this  case  two  different  lengths 
of  rough  sea  surface  are  considered  instead  of  flat  surfaces.  The  wind  speed  used  for  this 
realization  is  5  m/s  (RMS  wave  height  =  0.135  m).  The  rough  sea  surface  is  the  same  for 
both  cases  over  the  340  m  closest  to  the  target.  The  two  results  match  except  for  very  low 
angles.  The  results  of  this  figure  and  the  last  show  why  a  large  surface  must  be  considered 
when  low- grazing  angles  are  analyzed. 

Figure  2.23  shows  the  CPU  time  per  iteration  required  by  both  methods  (GFB  and  GFB 
with  Spectral  Acceleration)  as  a  function  of  the  number  of  unknowns  for  the  sea  surface. 


41 


Figure  2.22:  Bistatic  RCS  for  two  different  lengths  of  rough  sea  surface.  Frequency  =  300  MHz, 
elevation  angle  =  5°,  wind  speed  =  5  m/s. 

These  results  agree  with  the  computational  cost  estimate  shown  in  the  previous  section;  the 
GFB  method  presents  a  computational  cost  per  iteration  of  0(N2),  while  with  the  spectral 
acceleration  this  cost  is  reduced  to  O(N). 

2.4.5  Conclusions 

By  including  this  Spectral  Acceleration  in  the  Generalized  Forward-Backward  method,  it 
is  possible  to  analyze  the  scattering  from  targets  on  very  large  surfaces  with  a  low  compu¬ 
tational  cost  and  without  lost  of  accuracy.  This  makes  this  method  suitable  for  analyzing 
low-grazing  angle  scattering  problems  where  a  large  surface  must  be  considered  to  solve  the 
problem. 
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Chapter  3 

Plane  Wave  Spectral  Properties  of  the 
Sea  Scattered  Field 

3 . 1  Introduction 

In  this  chapter  the  plane  wave  spectrum  of  the  sea  scattered  field  due  to  an  incident  plane 
wave  is  investigated.  An  understanding  of  the  plane  wave  spectrum  is  important  for  gaining 
insight  into  a  ray-optical  representation  of  the  sea  scattered  field.  This  would  allow  a  target  to 
be  placed  on  the  sea  surface  and  its  RCS  computed  using  the  ray-optical  representation  as  the 
incident  field.  RCS  prediction  codes  which  use  high-frequency  asymptotic  methods  generally 
must  assume  a  ray-optical  incident  field.  In  the  following  the  polarization  is  horizontal  and 
normal  to  the  2D  geometry,  and  the  surface  is  assumed  to  be  PEC. 

The  plane  wave  spectrum  is  defined  from  the  Fourier  transform  of  the  fields  over  a  vertical 
cross-section  above  the  sea  surface,  as  shown  in  Figure  3.1.  The  total  forward  propagating 
electric  field  in  the  vicinity  of  the  target  is  given  approximately  in  terms  of  its  plane  wave 
spectrum  as 

/OO 

Fy(kz)  e-j{kxX+k‘z)dkz  (3.1) 

-OO 

where  the  Fourier  spectrum  is  given  by 

F,(k,)  =  ±J*Ey(0,z)e>k-’dz  (3.2) 

and  h  is  the  maximum  height  over  which  the  expansion  is  valid.  (3.1)  is  exact  over  the 
vertical  cross-section  at  x  =  0  and  becomes  more  approximate  as  x  is  varied  in  the  target 
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Figure  3.1:  Incident  and  sea  scattered  fields  transformed  into  a  spectrum  of  plane  waves  over 
a  vertical  cross-section  above  the  sea  surface. 

vicinity.  We  are  usually  only  interested  in  propagating  plane  waves  which  are  defined  in 
terms  of  a  propagation  angle  (f>  defined  as  the  angle  the  plane  wave  makes  with  the  x-axis. 
Using  the  transformation  kz  =  k  sin  <j>  and  integrating  only  over  the  propagating  spectrum, 
(3.1)  becomes 


Ey(x,z)  *  k  J*'*  Fy (0)  cos  4>  e~j(kx cos *+kz sin (3.3) 

where 

FM  =  J-  /* Et(0,z)  (3.4) 

Ztt  Jo 

This  approximation  is  expected  to  be  very  accurate  for  large  h  and  low  angles  of  incidence.  In 
the  results  which  follow  the  sea  scattered  fields  are  computed  numerically  using  the  forward- 
backward  method  described  in  Chapter  2.  The  fields  are  then  transformed  using  the  above 
expressions. 


3.2  Numerical  Results 

Figure  3.2  shows  the  magnitude  of  the  plane  wave  spectrum  for  a  5°  incidence  angle  and 
a  wind  speed  of  5  m/s  (11  mph)  at  1  GHz.  At  this  wind  speed  the  standard  deviation  of 
the  surface  height  profile  a  (RMS  wave  height)  is  0.135  m,  corresponding  to  Sea  State  2 
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as  defined  by  the  World  Meteorological  Organization  (WMO)  [2].  In  terms  of  wavelength, 
aj\  =  0.45.  In  general,  the  RMS  wave  height  in  meters  is  empirically  related  to  the  wind 
speed  u  in  meters/second  by  [2] 

a  =  0.0054m2.  (3.5) 

The  actual  surface  is  shown  at  the  top  of  the  figure.  The  “One  realization”  plot  shows  the 
total  plane  wave  spectrum  for  this  surface  computed  from  (3.4),  including  the  incident  plane 
wave  which  corresponds  to  the  impulse  at  a  scattering  angle  of  0  =  —5°.  The  height  h  over 
which  the  fields  are  transformed  for  these  results  is  100  m.  The  sea  scattered  spectrum  is 
visible  for  positive  scattering  angles,  and  shows  an  impulse  at  (j>  =  5°  which  corresponds  to 
the  coherent  reflection.  The  coherent  component  only  contributes  to  the  sea  scattered  field 
in  the  specular  (direct  reflection)  direction,  whereas  the  incoherent  component  contributes  in 
all  directions  [1].  The  coherent  reflection  is  expected  to  always  be  about  the  same  for  a  given 
incidence  direction  and  wind  speed,  independent  of  time.  The  use  of  the  term  “coherent” 
in  this  report  refers  to  a  quantity  that,  by  definition,  does  not  change  with  time.  This  is 
shown  in  the  “Coherent  part”  of  Figure  3.2,  which  is  the  total  spectrum  averaged  over  40 
different  randomly  generated  sea  surfaces  at  a  wind  speed  of  5  m/s.  The  incoherent  spectrum 
averages  out  to  zero  while  the  coherent  component  stays  relatively  constant. 

The  coherent  plane  wave  component  of  the  scattered  field  may  be  obtained  approximately 
by  the  physical  optics  (PO)  technique  used  in  [1]  as 

Ecy  «  El  T(^) 

r(^)  =  R(f)e-2(fe<TSin^i)2  (3.6) 

where  R{4>1)  is  the  plane  wave  reflection  coefficient  for  a  planar  interface  of  the  same  material 
as  the  sea  surface  and  4 1)1  is  the  elevation  angle  of  the  incidence  plane  wave.  For  the  case  of 
Figure  3.2,  |r|=0.89  computed  from  (3.6);  this  appears  to  agree  with  the  numerical  results 
because  the  coherent  sea  scattered  impulse  at  5°  is  about  0.89  of  the  incident  impulse  at 
-5°. 
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Figure  3.2:  Plane  wave  spectrum  due  to  a  1  GHz  plane  wa’ 
surface  with  wind  speed  5  m/s.  Coherent  and  incoherent  parti 


As  the  coherent  plot  of  Figure  3.2  shows,  the  remainder  of  the  total  spectrum,  the 
incoherent  spectrum,  averages  out  to  zero  over  many  surface  realizations  indicating  that 
it  has  a  uniformly  distributed  random  phase.  The  “Incoherent  part”  of  the  figure  plots 
the  RMS  power  of  the  incoherent  spectrum  averaged  over  40  different  randomly  generated 
surfaces  (and  plotted  on  a  smaller  scale).  For  this  result  the  coherent  spectrum  has  been 
subtracted  out.  The  incoherent  power  spectrum  appears  to  have  a  well-behaved  probability 
distribution  as  a  function  of  angle.  An  analytic  expression  for  this  distribution  is  being 
investigated  based  on  the  derivation  in  [1]. 

Figure  3.3  shows  the  spectra  for  the  same  plane  wave  incident  on  a  sea  surface  with  a 
wind  speed  of  10  m/s  (23  mph)  and  a  corresponding  RMS  wave  height  a= 0.54  m  (WMO 
Sea  State  4),  and  a/ A  =  1.8.  Clearly,  the  coherent  component  of  the  scattered  field  is 
much  weaker  than  before  (the  incident  plane  wave  component  stays  the  same).  For  this  case 
|pj— 0.14  from  (3.6),  which  is  in  fair  agreement  with  the  numerical  results.  Figure  3.4  shows 
the  spectra  for  a  sea  surface  with  a  wind  speed  of  15  m/s  (34  mph),  o= 1.22  m  (WMO  Sea 
State  5),  cr/A  =  4.1,  and  |r|=4.9  x  1(T5  from  (3.6).  For  this  case  the  coherent  sea  scattered 

component  is  not  discernible. 

As  Equation  (3.6)  indicates,  the  coherent  component  is  larger  for  lower  angles  of  incidence. 
Figure  3.5  shows  the  spectra  for  a  1  GHz  plane  wave  incident  at  1.5°  elevation  over  a 
sea  surface  at  15  m/s  wind  speed,  which  corresponds  to  |r|=0.41.  The  coherent  scattered 
component  is  again  visible  in  the  sea  scattered  spectrum.  This  indicates  that  the  coherent 
component  may  dominate  for  low-grazing  angle  RCS  problems.  This  will  be  investigated 

further  in  the  next  two  chapters. 
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Figure  3.3:  Plane  wave  spectrum  due  to  a  1  GHz  plane  wave  incident  at  5°  above  a  sea 
surface  with  wind  speed  10  m/s.  Coherent  and  incoherent  parts  are  obtained  by  sampling 
40  surface  realizations. 
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Wind  speed  =  1 5  m/s 
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0  (DEG) 


Figure  3.4:  Plane  wave  spectrum  due  to  a  1  GHz  plane  wave  incident  at  5°  above  a  sea 
surface  with  wind  speed  15  m/s.  Coherent  and  incoherent  parts  are  obtained  by  sampling 
40  surface  realizations. 
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Next,  the  plane  wave  reconstruction  of  the  fields  over  the  vertical  cross-section  above  the 
sea  surface  will  be  investigated.  Figure  3.6  plots  the  total  fields  and  the  fields  reconstructed 
using  only  the  coherent  components,  as  a  function  of  the  vertical  height  above  the  surface 
for  a  1  GHz  plane  wave  incident  at  5°.  Results  for  wind  speeds  of  5  m/s  and  10  m/s  are 
shown.  The  coherent  components  reconstruct  well  the  total  fields  for  the  5  m/s  case,  but 
not  for  the  10  m/s  case.  This  is  expected  from  a  comparison  of  Figures  3.2  and  3.3.  At 
5  m/s  the  coherent  scattered  component  is  much  greater  than  the  incoherent  spectrum,  so 
most  of  the  scattered  power  is  coherent.  But  at  10  m/s  the  coherent  scattered  component  is 
comparable  to  the  incoherent  spectrum,  so  more  power  is  distributed  in  the  incoherent  field. 

Since  it  is  clear  that  in  some  cases  the  coherent  field  is  insufficient  for  representing  the 
total  field,  the  following  discrete  plane  wave  expansion  is  considered: 

M 

Ey(x,z)  »  £  (3.7) 

771—1 

The  discrete  plane  wave  directions  (j)m  and  the  coefficients  Am  are  chosen  to  give  the  best 
reconstruction  for  a  given  number  of  plane  waves  M .  The  number  of  plane  waves  should  be 
as  small  as  possible  since  the  goal  of  the  plane  wave  expansion  is  to  provide  a  ray-optical 
incident  field  for  RCS  codes,  and  efficiency  is  a  primary  concern.  It  would  be  inefficient,  for 
example,  to  sample  the  continuous  spectrum  of  (3.3)  at  the  Nyquist  rate  and  replace  the 
integral  with  a  summation.  This  is  because  the  height  h  over  which  the  expansion  is  needed 
to  illuminate  the  target  could  be  quite  large  in  terms  of  wavelength,  so  a  very  large  number 
of  plane  waves  would  be  needed.  Instead,  it  is  possible  to  extract  only  the  largest  plane  wave 
components  from  a  given  sea  scattered  spectrum  using  a  simple  feature  extraction  algorithm. 
Figure  3.7  shows  the  reconstruction  of  the  10  m/s  field  using  the  50  strongest,  plane  wave 
components,  instead  of  just  the  coherent  components  as  in  Figure  3.6.  The  reconstruction 
agrees  much  better  with  the  total  field.  This  discrete  plane  wave  expansion  will  be  used  as 
the  incident  field  for  computing  the  RCS  of  a  target  on  the  surface  in  Chapter  4. 
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Wind  speed  =  5  m/s 


Wind  speed  =  10  m/s 


- Total  field  - Coherent  components  only  (including  incident  plane  wave) 


Figure  3.6:  Coherent  plane  wave  reconstruction  of  the  total  fields  above  the  sea  surface. 
Top:  wind  speed  =  5  m/s.  Bottom:  wind  speed  =  10  m/s.  Incident  elevation  angle  =  5°, 
frequency  =  1  GHz. 
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height  (m) 

- Total  field  - 50  strongest  plane  waves  (plus  incident  plane  wave) 

Figure  3.7:  Reconstruction  of  the  total  fields  above  the  sea  surface  using  the  50  strongest 
plane  wave  components.  Wind  speed  =  10  m/s.  Incident  elevation  angle  =  5°,  frequency  = 
1  GHz. 

3.3  Conclusions 

It  has  been  shown  that  the  plane  wave  spectrum  of  the  sea  scattered  field  due  to  an  incident 
plane  wave  may  be  decomposed  into  coherent  and  incoherent  components.  The  coherent 
component  is  well  predicted  by  the  PO  formula  of  (3.6),  and  may  be  stronger  or  weaker 
than  the  incoherent  spectrum  depending  on  the  incidence  angle,  frequency,  and  wave  height. 
As  the  incidence  angle  approaches  grazing,  the  coherent  component  strength  approaches  the 
strength  of  the  geometrical  optics  reflected  field  for  a  planar  interface,  so  by  conservation 
of  power  the  incoherent  component  goes  to  zero.  For  shipboard  surface-to-surface  radar 
systems,  the  incidence  angle  is  almost  always  very  near  grazing  because  the  transmit  antenna 
height  is  much,  much  smaller  than  the  distance  to  the  target.  Even  for  air-to-surface  radars 
the  incidence  angle  is  probably  going  to  be  very  close  to  grazing  for  most  problems  of 
interest.  Therefore,  in  most  cases  the  coherent  component  is  probably  going  to  be  sufficient 
for  representing  the  sea  scattered  field  in  the  vicinity  of  the  target. 

For  cases  when  the  coherent  component  alone  does  not  accurately  represent  the  sea 
scattered  field,  a  discrete  plane  wave  expansion  may  used  as  in  (3.7).  A  relatively  small 
number  of  plane  waves  may  be  chosen  if  the  strongest  spectral  components  can  be  found  from 


a  given  spectrum.  The  possibility  of  using  discrete  plane  waves  with  randomly  generated 
amplitudes  and  phases  is  being  investigated,  and  is  based  on  ascertaining  the  probability 
distribution  function  of  the  sea  scattered  spectrum  for  a  given  wind  speed.  The  effect 
of  using  the  plane  wave  expansion  to  illuminate  a  target  on  a  surface  for  the  purpose  of 
computing  its  RCS  is  investigated  in  the  next  chapter. 
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Chapter  4 


Ship  Radar  Cross  Section 
Computation  Using  Ray-Optical 
Methods 

4.1  Introduction 

In  this  chapter  the  RCS  of  a  2D  target  on  a  sea  surface  is  computed  using  ray-optical  methods 
that  decouple  the  ship  scattering  problem  from  the  sea  surface  scattering.  First,  the  total 
incident  field  in  the  vicinity  of  the  target  (i.e.,  incident  plane  wave  plus  sea  scattered  field), 
but  with  the  target  absent,  is  expressed  in  terms  of  a  discrete  set  of  plane  waves  as  in  the  last 
chapter.  Next,  the  target  is  placed  on  an  infinite  planar  interface  and  illuminated  with  the  set 
of  plane  waves  to  find  the  scattered  fields  over  the  surface  of  the  target.  Lastly,  the  reciprocity 
theorem  is  used  to  find  the  backscattered  field  by  reacting  the  incident  field  with  the  scattered 
field  over  the  surface  of  the  target.  The  reciprocity  formulation  is  derived  first  in  Section  4.2 
because  it  defines  the  field  quantities  needed  to  compute  the  RCS.  Then  in  Section  4.3  the 
approximately  equivalent  ray-optical  scattering  problem  is  presented.  Numerical  results  are 
presented  which  compare  the  ray-optical  results  with  rigorous  numerical  results  generated 
using  the  spectrally  accelerated  GFB  method.  Conclusions  follow  at  the  end  of  the  chapter. 
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4.2  Reciprocity  Formulation  for  the  Backscattered  Field 

4.2.1  Basic  Derivation 

The  following  derivation  is  for  a  three-dimensional  geometry  except  where  noted. 

Figure  4.1  shows  three  different  field  configurations  over  the  same  sea  surface.  The  top 
figure  is  the  geometry  of  interest,  which  is  a  scattering  target  on  the  sea  surface  illuminated 
by  an  incident  field  generated  by  the  antenna  on  the  left,  and  the  desired  quantity  is  the 
scattered  field  ( ES,HS )  due  to  the  presence  of  the  target.  The  incident  field  (E\Hl)  is 
defined  as  the  field  over  the  sea  surface  in  the  absence  of  the  ship,  and  is  illustrated  in  the 
middle  figure.  The  total  field  is  therefore, 


(E,H)=(Ei,Hi)+(Es,Hs) 


Incident  Fields 
w/o  ship 


Test  Fields 
w/o  ship 


Figure  4.1:  Sea  surface  and  scattering  target  configurations  for  the  reciprocity  formulation. 
Top:  incident  and  scattered  fields  with  target  present.  Middle:  Incident  fields  with  target 
absent.  Bottom:  Fields  of  a  test  source  with  target  absent. 


(E,H)  =  {E^H^  +  iE^H8). 


(4.1) 
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It  is  of  interest  to  compute  the  scattered  electric  field  at  some  point  Pt.  To  invoke  reciprocity, 
we  place  an  electric  point  dipole  test  source  dpt  at  Pt  and  define  the  fields  (Et,  Ht )  as  the 
fields  of  this  test  source  in  the  presence  of  the  sea  surface,  but  in  the  absence  of  the  target, 
as  illustrated  in  the  bottom  figure  of  Figure  4.1.  Applying  the  reciprocity  theorem  [4]  to  the 
fields  of  the  test  source  and  the  scattered  fields  within  the  volume  V  yields 

dpt  ■  Es(Pt)  =  f  (Es  xHt-EtxHs)-ndS  (4.2) 

JSsea+Soo+ST  V  ' 

where  h  points  into  the  volume  V.  The  volume  V  contains  the  point  Pt  and  is  bounded  by 
the  target  surface  St,  the  sea  surface  Ssea  outside  of  the  target,  and  the  surface  at  infinity 
Soo.  The  fields  ( Es ,  Hs)  and  (Et,  Ht )  satisfy  the  same  boundary  conditions  on  Ssea  and  S 
so  the  surface  integral  vanishes  on  those  surfaces  leaving 

dpt  ■  Es(Pt)  =  f  (Es  x  Ht-EtxHs)-hdS  (4.3) 

J  St 

which  is  the  form  of  interest.  This  equation  states  that  the  scattered  electric  field  at  Pt  can 
be  found  by  integrating  the  “reaction”  of  the  scattered  fields  with  the  test  fields  over  the 
target  surface.  (The  reciprocity  integral  is  sometimes  called  the  reaction  integral.)  This  is 
computationally  significant  because  the  target  scattered  fields  only  need  to  be  found  on  the 
surface  of  the  target,  and  do  not  need  to  be  tracked  from  the  target  back  to  the  observer. 
The  expression  of  (4.3)  is  an  exact  formulation  derived  directly  from  Maxwell’s  equations 
and  from  the  boundary  conditions  on  the  sea  surface  and  at  infinity.  An  approximation  will 
be  used  later  to  compute  the  scattered  fields  ( ES,HS )  on  the  surface  of  the  target. 


4.2.2  Simplification  for  Plane  Wave  Backscattering 

For  backscattering  problems  the  formulation  is  even  simpler  because  the  incident  source  may 
be  used  as  the  test  source,  so  (Et,Ht)  —  ( El,Hl ).  A  plane  wave  incident  field  is  generated 
by  placing  a  dipole  source  at  infinity,  as  shown  in  Figure  4.2.  The  dipole  p  is  assumed  to 
have  unit  strength  and  is  placed  at  f,  r  — >  oo  with  respect  to  the  reference  origin  O.  The 
plane  wave  incident  electric  field  at  a  O  as  shown  in  the  figure  is  then  given  by 

El0(0)  =  Urn  -p  jkZ°e  (4.4) 

u  v  7  r— >oo  47T  T 
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(E,H)=(Ei,Hi)+(Es,Hs) 


(E1,!!1)  are  in 
presence  of  sea 
surface 


Figure  4.2:  Plane  wave  backscattering  configuration  for  the  reciprocity  formulation. 


The  reciprocity  equation  then  becomes 

p  ■  Es(f)  =  lim  -^2  e~jkr  [  (Es  x  it*  -  &  x  Hs)  ■  h  dS  (4.5) 

where  (El,Hl)  are  the  fields  due  to  the  incident  plane  wave,  phase-referenced  to  point  O, 
in  the  presence  of  the  sea  surface.  This  equation  may  be  used  directly  for  backscatter  RCS 
computations.  It  is  only  necessary  to  find  the  scattered  fields  ( Es ,  Hs)  over  the  target  surface 
due  to  the  incident  fields  ( El,Hl ).  For  2D  scattering  problems,  a  line  source  replaces  the 
point  source  and  the  factor  is  replaced  with  Z0^j ^ 


4.2.3  Equivalent  Current  Form  for  a  PEC  Target 


Equivalent  surface  currents  are  often  used  in  RCS  computation  methods,  such  as  in  physical 
optics  (PO)  [3]  and  the  method  of  moments  [7].  For  a  PEC  target,  the  equivalent  magnetic 
surface  current  is  zero  and  the  electric  current  is  given  by 


Js  =  n  x  H 


(4.6) 


where  H  is  the  total  field.  In  the  integrand  of  (4.5),  for  a  PEC  target  it  may  be  shown  that 


(£s  xHi-EixHs)-h  = 


-E{  xHi-Ei  x(H 


=  -EixH-n 
=  El  •  Js 


(4.7) 


so  (4.5)  becomes 

p  •  Es(f)  =  lim  e~jkr  [  E{  ■  Js  dS. 

r->oo  4tt  r  JsT 


(4.8) 


59 


The  problem  is  reduced  to  finding  Js  over  the  target  surface.  It  must  be  noted  that  Js 
is  found  in  the  presence  of  the  sea  surface,  so  the  usual  equivalent  current  methods  must 
be  modified.  For  example,  the  electric  field  integral  equation  (EFIE)  [7]  for  this  problem  is 
given  by 

-fix  E*  =  fix  [  Js-  Gsea  dS  (4.9) 

J  St 

where  Gsea  is  the  dyadic  Green’s  function  for  the  space  above  the  sea  surface  instead  of  the 
usual  free  space  dyadic  Green’s  function  G0. 

4.2.4  Local  Planar  Approximation  for  Sea  Surface 

Since  the  Green’s  function  for  a  rough  surface  is  not  easy  to  find,  an  approximation  may  be 
used  which  assumes  the  surface  is  locally  flat  only  for  the  purposes  of  finding  Js  in  (4.9). 
The  incident  field  E{  is  still  in  the  presence  of  a  rough  surface.  With  this  approximation, 
the  sea  surface  dyadic  Green’s  function  in  (4.9)  is  replaced  with: 

Gsea  ~  G0  +  Gi  (4-10) 


where  Gj  is  the  image  dyadic  Green’s  function  as  shown  in  Figure  4.3.  Using  image  theory 
[4],  the  image  of  the  incremental  current  element  generates  the  reflected  field.  It  is  noted 
that  Gi  must  take  into  account  the  reflection  coefficient  from  the  planar  interface  if  it  is 
not  PEC.  The  image  approach  is  not  exact  for  non-PEC  interfaces,  but  should  be  a  good 
approximation  if  the  interface  is  of  a  highly  impenetrable  material  such  as  sea  water.  The 
reciprocity  integral  for  plane  wave  incidence  (4.9)  becomes 

—h  x  Ei  =  fix  [  (Js-Go  +  Js-Wj)  dS.  (4.11) 

J  St 

4.3  Equivalent  Ray-Optical  Scattering  Problem 

In  this  section  the  discrete  plane  wave  expansion  of  (3.7)  in  Chapter  3  is  used  to  define  the 
incident  field  in  the  local  planar  approximation  for  the  EFIE  (4.11)  for  finding  the  equivalent 
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Source 


Figure  4.3:  Local  planar  approximation  for  the  sea  surface  and  image  theory  used  to  find 
the  scattered  fields  in  the  vicinity  of  the  target. 


currents  over  the  target  surface.  Referring  to  Figure  4.4,  the  incident  fields  at  a  point  r'  on 
the  target  surface  are  given  by 


E‘(f) 

Hi(f ) 


M 

£  pm  Ame~jkm'r 


m= 0 


M  y 

^  ^  km  ^  Pm  f  ry  AmC  j 
m= 0 


kZn 


(4.12) 


where  for  the  mth  plane  wave,  pm  is  the  unit  electric  polarization  vector,  km  is  the  propagation 
vector  (|A;m|  =  k),  and  Am  is  the  bistatic  sea  scattering  coefficient  in  the  direction  of  km. 
The  m  =  0  term  is  the  original  incident  plane  wave,  and  the  m  >  0  terms  are  the  sea 
scattered  plane  waves.  For  the  case  when  only  the  coherent  components  are  used,  M  —  1 
and  k\  =  k0  —  2 (z  •  kQ)z  (where  z  points  straight  up),  and  Ai  =  T(k0)  from  Equation  (3.6). 

Using  the  above  ray-optical  incident  field  expansion,  the  EFIE  of  (4.11)  is  solved  numeri¬ 
cally  for  Js  over  the  ship  surface  using  the  MoM  [7]  or  approximately  using  a  high-frequency 
method  such  as  physical  optics  [3].  It  is  also  possible  to  use  (4.5)  directly  to  find  the  scat¬ 
tered  fields  using  high-frequency  methods,  without  the  intermediate  step  of  first  solving  for 
equivalent  currents.  In  the  numerical  results  of  the  next  section  the  MoM  is  used.  It  is  noted 
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(E^)  ° 


Figure  4.4:  Target  illuminated  by  a  discrete  set  of  plane  waves. 

again  that  (4.11)  is  solved  in  the  presence  of  a  locally  flat  surface  using  the  image  Green’s 
function.  This  allows  the  target  scattering  problem  to  be  solved  approximately,  independent 
of  the  exact  rough  sea  surface. 

As  an  example  of  a  high-frequency  solution  using  physical  optics  [3],  the  surface  currents 
may  be  found  approximately  by 

Js(f')  «  2 nxH'GO(f')  (4.13) 


where  HzGO(r')  is  the  geometrical  optics  (GO)  component  of  the  incident  field.  In  terms  of 
the  plane  wave  expansion  of  (4.12),  the  PO  currents  are  given  by 


Ur') 


£ 


m=0 


ft  X  (krn  X  Pm)  kZ0^m^  ^  ^ 
0 


lit  region  of  mth  plane  wave, 
shadow  region  of  mth  plane  wave. 

(4.14) 


In  the  above  the  “lit”  region  is  defined  as  points  f'  on  the  surface  of  the  target  where  the 
GO  field  is  non-zero.  As  with  the  EFIE,  the  GO  field  must  also  include  the  locally  planar 
surface  during  the  ray  tracing  process,  as  illustrated  in  Figure  4.5.  It  is  possible  that  a  GO 
ray  may  reflect  from  the  target  and  hit  the  sea  surface,  then  hit  the  target  again  as  shown. 
The  second  hit  on  the  target  would  also  contribute  to  3S  in  the  PO  approximation.  This 
contribution  is  important  to  include  because  the  ship  hull  and  its  image  could  form  a  large 
corner  reflector,  which  would  likely  be  the  dominant  backscatter  mechanism  for  the  target. 
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Figure  4.5:  GO  ray  from  the  mth  incident  plane  wave  which  reflects  from  the  target,  hits  the 
sea  surface,  and  then  hits  the  target  again. 

4.4  Numerical  Results 


In  the  numerical  results  of  this  section  the  general  3D  formulations  of  the  previous  section  are 
converted  to  2D,  and  the  electric  field  is  horizontally  polarized  (i.e.,  pm  =  y  and  pm  •  km  =  0) 
as  in  Chapters  2  and  3.  The  equivalent  surface  currents  on  the  target  are  found  using  (4.11) 
with  (4.12)  representing  the  incident  fields,  and  the  backscattered  fields  are  found  using 
(4.8).  The  2D  RCS,  or  more  specifically,  the  2D  radar  echo  width,  is  defined  by 

£J(f)f.  (4.15) 


RCS2£>  =  lim  27rr 


Using  the  2D  form  of  (4.8)  for  horizontal  polarization  this  becomes 

kZl  '  '  12 


RCS 


2D 


=  zzi  f 

4  J St 


Ely  Js  dS 


(4.16) 


where  Js  is  the  scalar  portion  of  Js  which  is  y-directed. 

The  frequency  is  1  GHz  in  all  the  results  of  this  section.  The  reference  RCS  results 
predicted  by  the  GFB  method  are  for  the  total  backscatter  from  the  ship  and  sea  surface 
together.  However,  as  shown  in  the  results  of  Chapter  2,  the  backscatter  from  the  sea  surface 
is  generally  more  than  30  dB  below  the  backscatter  from  the  ship. 

To  establish  an  accuracy  bound  on  the  GFB  reference  solution,  Figure  4.6  plots  the  2D 
RCS  of  the  13  m  high  ship  target  on  a  402.4  m  flat  sea  surface  computed  using  the  spectrally 
accelerated  GFB  method.  The  2D  ship  geometry  is  shown  in  the  figure  above  the  plot;  the 
hull  makes  an  angle  of  80°  with  the  horizontal  axis.  The  plane  wave  is  incident  from  the  left. 
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Elevation  angle  (degrees) 


Figure  4.6:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  flat  surface  computed  using  GFB  and 
image  theory.  Frequency  =  1  GHz. 


Also  plotted  is  a  reference  solution  computed  using  exact  image  theory  for  the  ship  on  an 
infinite  flat  surface.  It  is  noted  that  for  a  flat  sea  surface  the  sea  scattered  incident  field  is 
entirely  coherent.  Notice  the  large  peak  around  10°  which  is  due  to  the  sea  reflected  plane 
wave  hitting  the  ship  hull  broadside.  There  is  excellent  agreement  everywhere  except  below 
about  3°.  This  is  due  to  the  fact  that  the  GFB  uses  a  finite  flat  surface  (with  R-cards  on  the 
ends),  whereas  the  image  theory  result  assumes  an  infinite  surface.  The  GFB  solution  could 
be  improved  for  low  angles  by  using  a  larger  sea  surface,  but  at  the  expense  of  increased 
CPU  time.  As  an  example  of  the  CPU  time  required,  the  spectrally  accelerated  GFB  result 
took  2  hours  to  compute  on  a  450  MHz  Pentium  II  processor,  and  the  image  theory  result 
about  2  minutes.  The  patterns  are  sampled  at  0.2°  increments.  These  times  are  typical  of 
all  the  results  of  this  chapter — 2  hours  for  GFB  RCS  patterns  and  2  minutes  for  coherent 
RCS  patterns. 

Figure  4.7  plots  the  2D  RCS  of  the  13  m  high  ship  target  on  a  402.4  m  sea  surface  with 
wind  speed  5  m/s  (RMS  wave  height  a  =  0.14  m,  a/\  =  0.47).  The  ship  geometry  and 
the  last  100  m  of  sea  surface  are  shown  above  the  plot.  The  solid  line  is  the  GFB  reference 
solution,  and  the  dashed  line  is  the  RCS  computed  using  only  the  coherent  components  of 
the  incident  field  (in  the  presence  of  the  sea  surface)  and  uses  the  local  planar  approximation 
of  (4.11)  solved  using  the  MoM.  The  coherent  component  of  the  sea  scattered  incident  field 
is  found  from  the  formula  of  (3.6).  The  agreement  is  good  for  low  angles  (except  for  the 
discrepancies  at  very  low  angles  due  to  the  finite  surface),  but  the  coherent  result  generally 
under-predicts  the  RCS  at  higher  angles.  This  is  expected  from  (3.6)  since  the  coherent 
component  becomes  weaker  as  the  elevation  angle  increases. 

Figure  4.8  shows  the  same  RCS  result,  but  the  dashed  line  is  computed  using  the  50 
strongest  plane  wave  components  to  represent  the  sea  scattered  incident  field  (in  addition 
to  the  incident  plane  wave  k0).  The  plane  wave  amplitudes  are  found  by  extracting  the 
strongest  spectral  components  of  the  sea  scattered  field  with  the  ship  absent,  which  was 
computed  using  the  FB  method.  The  overall  agreement  with  the  reference  GFB  solution 
does  not  improve  much  compared  with  the  coherent  result,  although  the  envelope  of  the 
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Figure  4.7:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  surface  computed  using  coherent 
components  of  sea  scattered  incident  field.  Wind  speed  =  5  m/s,  frequency  =  1  GHz. 


66 


patterns  appears  to  agree  better  for  higher  angles.  The  RCS  is  generally  not  under-predicted 
in  this  plot,  except  for  a  few  isolated  angles.  This  is  important  because  the  purpose  of  the 
approximate  ray-optical  methods  is  to  predict  the  peak  levels  and  range  of  variation  of  the 
RCS  pattern. 

Figure  4.9  shows  the  RCS  computed  using  the  numerically  exact  incident  field  obtained 
from  the  FB  method.  The  results  are  very  similar  to  the  50  plane  wave  incident  field 
representation,  indicating  again  that  the  plane  waves  adequately  reconstruct  the  sea  scattered 
incident  field.  The  differences  with  the  GFB  reference  solution  are  due  entirely  to  the  local 
planar  approximation  used  to  find  the  scattered  fields  over  the  surface  of  the  target,  which 
seems  to  mostly  effect  the  details  of  the  pattern  and  not  the  overall  envelope. 

Figure  4.10  shows  the  RCS  of  the  13  m  ship  computed  using  the  GFB  method  for  a  flat 
surface  and  for  Surface  #1  (wind  speed  =  5  m/s).  Also  plotted  is  the  RCS  computed  using 
only  the  coherent  component  of  the  incident  field.  This  figure  shows  that  the  envelope  of 
the  RCS  pattern  for  the  rough  surface  is  generally  bounded  above  by  the  envelope  of  the  flat 
surface  pattern,  and  below  by  the  envelope  of  the  coherent  pattern.  Also  notice  that  the  two 
GFB  results  agree  very  well  for  low  angles  because  they  both  include  the  effects  of  using  a 
finite  sea  surface  for  the  computations. 

To  see  if  the  envelope  agreement  holds  up  for  a  different  randomly  generated  surface 
with  the  same  wind  speed,  Figure  4.11  shows  the  same  results  for  Surface  #2  which  is  shown 
above  the  plot.  The  two  GFB  results  agree  quite  well  for  this  case,  and  the  envelopes  of  the 
flat  and  coherent  results  provide  an  upper  and  lower  bound  for  the  rough  surface  RCS  plot. 
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Figure  4.9:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  surface  computed  using  the  numerically 
exact  sea  scattered  incident  field.  Wind  speed  =  5  m/s,  frequency  =  1  GHz. 


Elevation  angle  (degrees) 

Figure  4.10:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  surface  computed  using  GFB  for  a 
flat  surface  and  a  rough  surface.  Wind  speed  =  5  m/s,  frequency  =  1  GHz. 
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Figure  4.11:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  surface  computed  using  GFB  for  a 
flat  surface  and  a  rough  surface.  Wind  speed  =  5  m/s,  frequency  =  1  GHz. 


To  investigate  further  the  envelope  behaviors  of  the  RCS  patterns,  Figure  4.12  shows  the 
RCS  of  the  13  m  ship  on  a  surface  with  a  wind  speed  of  10  m/s  (RMS  wave  height  a  =  0.54  m, 
cr/A  =  1.8),  and  Figure  4.13  shows  the  RCS  for  a  wind  speed  of  15  m/s  (RMS  wave  height 
a  =  1.22  m,  a/X  =  4.1).  The  ship  geometry  (which  is  the  same)  and  the  last  200  m  of 
the  sea  surface  are  shown  above  the  plots.  Again,  the  flat  surface  and  coherent  results  are 
also  plotted.  The  envelope  of  the  rough  surface  RCS  patterns  generally  lie  between  the  flat 
surface  and  coherent  envelopes,  and  the  overall  level  of  the  flat  surface  pattern  is  in  good 
agreement  with  the  level  of  the  rough  surface  pattern.  The  patterns  generally  show  more 
variation  in  the  RCS  at  higher  wind  speeds,  especially  at  lower  elevation  angles. 

Figure  4.14  shows  the  same  results  as  Figure  4.13  but  for  a  different  randomly  generated 
surface  at  the  same  wind  speed  of  15  m/s.  Even  though  the  patterns  for  Surface  #4  and 
Surface  #5  are  significantly  different,  the  same  conclusions  about  the  envelope  levels  seem 
to  be  confirmed. 
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Figure  4.12:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  surface  computed  using  GFB  for  a 
flat  surface  and  a  rough  surface.  Wind  speed  =  10  m/s,  frequency  =  1  GHz. 


Figure  4.13:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  surface  computed  using  GFB  for  a 
flat  surface  and  a  rough  surface.  Wind  speed  =  15  m/s,  frequency  =  1  GHz. 
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Figure  4.14:  2D  RCS  of  a  13  m  high  ship  on  a  402.4  m  surface  computed  for  flat  case, 
coherent  case,  and  rough  surface  case.  Wind  speed  =  15  m/s,  frequency  =  1  GHz. 


All  of  the  previous  results  are  for  the  same  ship  with  a  hull  angle  of  80°  with  respect 
to  the  sea  surface.  It  may  be  noticed  that  the  peak  levels  and  the  variations  of  the  RCS 
patterns  for  this  target  at  different  wind  speeds  are  about  the  same,  even  for  zero  wind 
speed  (flat  surface).  Figure  4.15  plots  the  RCS  patterns  for  a  hull  angle  of  90°,  with  the 
rest  of  the  topside  ship  geometry  unchanged.  The  wind  speed  is  15  m/s,  and  the  sea  surface 
is  the  same  as  in  Figure  4.13  (Surface  #4).  The  90°  hull  angle  is  a  pathological  case  (and 
is  often  encountered)  because  the  side  of  the  ship  forms  a  large  corner  reflector  with  a  flat 
sea  surface,  so  the  RCS  is  expected  to  be  very  high.  The  GFB  Flat  Surface  result  confirms 
this.  Unlike  the  previous  target  with  an  80°  hull  angle,  the  rough  surface  greatly  reduces 
the  peak  RCS  levels  for  this  target.  The  rough  surface  result  is  bounded  by  the  flat  surface 
and  coherent  results  as  before,  although  the  space  between  them  is  much  larger. 

Lastly,  it  is  of  interest  to  investigate  the  use  of  a  high-frequency  asymptotic  technique  for 
computing  the  RCS  of  a  simple  target  geometry.  Figure  4.16  plots  the  RCS  of  a  flat-top  ship 
on  a  flat  surface  (shown  above  the  plot),  computed  using  the  EFIE  from  (4.11)  and  physical 
optics  from  (4.14).  It  is  noted  that  for  this  flat  surface  case  only  the  incident  and  direct 
reflected  plane  waves  are  used.  The  ship  is  9  m  high  and  the  hull  makes  an  angle  of  80°  with 
the  surface.  The  PO  solution  is  first  order,  so  multi-bounce  effects  are  not  included.  The 
agreement  between  the  numerically  exact  EFIE  solution  and  the  approximate  PO  solution 
is  very  good.  The  strong  peak  around  10°  is  due  to  the  sea-reflected  plane  wave  hitting  the 
hull  broadside. 
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Figure  4.15:  2D  RCS  of  a  13  m  high  ship  with  a  90°  hull  angle  on  a  402.4  m  surface  computed 
for  flat  case,  coherent  case,  and  rough  surface  case.  Wind  speed  =  15  m/s,  frequency  = 


Figure  4.16:  2D  RCS  of  a  9  m  high  ship  on  a  402.4  m  flat  surface  computed  with  the  EFIE 
and  physical  optics.  Frequency  =  1  GHz. 
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4.5  Conclusions 


The  numerical  results  show  that  computing  the  RCS  of  a  target  on  a  rough  surface  using 
only  the  incident  plane  wave  and  the  coherent  plane  wave  component  of  the  sea  scattered 
field  to  illuminate  the  target,  generally  under-predicts  the  actual  RCS  level.  Using  more 
plane  waves  to  characterize  the  sea  scattered  field  raises  the  RCS  level  to  where  it  should 
be,  but  the  details  of  the  pattern  are  not  recovered.  Of  course,  the  advantage  of  using  plane 
waves  to  characterize  the  sea  scattered  field  which  illuminates  the  ship  is  the  great  reduction 
in  CPU  time,  because  the  sea  surface  is  then  decoupled  from  the  ship  scattering  calculation. 
What  is  important  is  that  the  ray-optical  results  predict  the  peak  levels  of  the  RCS  patterns 
and  the  range  of  variation. 

The  numerical  results  also  show  that  the  RCS  of  a  target  on  a  rough  surface  is  generally 
bounded  above  by  the  RCS  of  the  same  target  on  a  flat  surface.  However,  as  in  the  case  with 
the  90°  hull  angle  (corner  reflector),  the  flat  surface  RCS  may  be  considerably  higher  than 
the  rough  surface  RCS.  And  in  the  next  chapter  it  will  be  seen  that  the  RCS  as  a  function  of 
time  on  a  moving  sea  surface  can  have  peaks  above  the  flat  surface  RCS  for  some  targets.  It 
is  shown  that  the  RCS  as  a  function  of  roll  angle  should  also  be  investigated  to  understand 
the  RCS  variations  of  a  target  on  a  realistic  sea  surface. 
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Chapter  5 

Statistical  Properties  of  the  Ship 
Radar  Cross  Section 

5.1  Introduction 

In  the  last  chapter  the  ray-optical  approaches  were  shown  to  be  useful  for  for  finding  the 
bounds  on  the  RCS  pattern  of  a  target  on  a  rough  sea  surface,  as  a  function  of  the  elevation 
angle  of  the  incident  plane  wave.  In  general,  the  pattern  is  bounded  below  by  the  coherent 
RCS  computation,  and  bounded  above  by  the  RCS  pattern  for  the  target  on  an  infinite  flat 
surface.  These  bounding  patterns  are  easy  to  find  because  they  only  require  a  simple  plane 
wave  scattering  analysis  combined  with  image  theory.  The  effects  of  differing  random  sea 
surfaces  and  wind  speeds  on  the  RCS  patterns  were  also  investigated. 

In  this  chapter  the  statistical  properties  of  the  RCS  of  a  ship-like  target  on  a  rough  sea 
surface  are  investigated.  For  these  results  the  elevation  angle  of  the  incident  plane  wave 
and  the  wind  speed  will  be  held  constant  while  other  parameters  are  varied.  In  particular, 
the  sea  surface  will  evolve  as  a  function  of  time  so  that  the  RCS  as  a  function  of  time  may 
be  viewed.  Also,  the  ship  roll  angle  will  be  varied.  Lastly,  different  target  shapes  will  be 
considered.  Conclusions  are  discussed  at  the  end  of  the  chapter. 
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5.2  RCS  as  a  Function  of  Time 


The  RCS  plotted  as  a  function  of  time  is  probably  the  most  realistic  visualization  of  how  the 
rough  sea  surface  affects  the  RCS.  In  the  results  which  follow,  an  initial  random  sea  surface 
is  allowed  to  evolve  with  time  according  to  a  deep-water  linear  hydrodynamic  dispersion 
relation.  The  RCS  is  computed  at  each  time  step  using  the  GFB  method.  The  ship  is 
allowed  to  float  up  and  down  with  the  waves  so  that  the  left  side  of  the  hull  is  always  the 
same  length  above  the  waterline  for  each  time  step,  but  the  ship  is  not  allowed  to  roll  as  a 
function  of  time.  (The  roll  of  a  ship  is  a  very  complex  function  of  the  moving  sea  surface, 
and  depends  on  the  ship  center  of  gravity,  inertial  moment,  etc.,  which  is  beyond  the  scope 
of  this  investigation.) 

Figure  5.1  shows  the  RCS  as  a  function  of  time  for  the  13  m  high  ship  with  a  90°  hull 
angle  of  Figure  4.15  of  the  last  chapter.  The  wind  speed  is  15  m/s,  and  at  time  zero  the 
sea  surface  is  exactly  the  same  as  Surface  #4  of  Figure  4.15.  The  elevation  angle  of  the 
incident  plane  wave  is  held  constant  at  5°.  Also  shown  is  the  RCS  level  for  the  ship  on  a  flat 
surface  and  the  coherent  RCS,  both  of  which  obviously  do  not  vary  with  time.  Somewhat 
surprisingly,  the  peaks  of  the  rough  surface  RCS  never  reach  the  flat  surface  level,  as  might 
be  expected  if  the  sea  surface  becomes  locally  flat  near  the  ship.  Figure  5.2  shows  the  actual 
evolution  of  Surface  #4  at  5  second  intervals.  The  vertical  scales  in  this  figure  are  stretched 
to  better  visualize  the  differences  in  the  surfaces.  At  time  zero  the  surface  appears  to  be 
relatively  flat  just  in  front  of  the  ship,  so  the  RCS  is  expected  to  be  high  because  of  the 
corner  reflector  effect.  The  RCS  pattern  shows  a  peak  at  zero,  but  it  is  considerably  lower 
than  the  flat  surface  level.  At  5  s  and  15  s  the  top  of  the  ship  reaches  its  highest  points  on 
top  of  the  waves,  but  the  RCS  does  not  show  peaks  there.  The  RCS  is  low  at  10  s  and  20  s, 
but  it  is  not  immediately  apparent  how  the  local  sea  surface  affects  this. 
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Figure  5.1:  2D  RCS  vs.  time  for  a  13  m  high  ship  on  a  402.4  m  surface.  Elevation  angle 
5°,  wind  speed  =  15  m/s,  frequency  =  1  GHz. 


Figure  5.2:  13  m  high  ship  on  a  402.4  m  time-evolving  sea  surface.  Wind  speed  =  15  m/s, 
frequency  =  1  GHz. 


Figure  5.3  shows  the  same  RCS  vs.  time  result,  with  additional  plots  for  wind  speeds  of 
5  m/s  and  10  m/s.  The  5  m/s  result  is  very  close  to  the  flat  surface  level,  and  the  10  m/s 
result  has  peaks  a  little  higher  than  the  15  m/s  result.  Furthermore,  the  5  m/s  result  does 
not  vary  much  from  its  mean,  while  the  10  m/s  and  15  m/s  results  vary  quite  significantly. 
The  15  m/s  result  is  the  most  rapidly  varying,  although  the  range  of  variation  is  similar  to 
the  10  m/s  result. 

Figure  5.4  shows  the  RCS  vs.  time  for  the  13  m  ship  with  a  hull  angle  of  80°  at  a  wind 
speed  of  15  m/s.  Also  shown  is  the  ship  RCS  on  a  flat  surface  and  the  coherent  RCS.  The 
levels  are  lower  than  the  90°  hull  angle  case,  but  the  peaks  of  the  time-varying  RCS  come 
closer  to  the  flat  surface  level.  The  time-varying  RCS  pattern  stays  generally  between  the 
flat  surface  and  coherent  levels. 

Figure  5.5  shows  the  backscattered  fields  as  a  function  of  time  for  a  directive  300  MHz 
source  antenna  which  is  150  m  from  the  ship,  as  shown  in  the  figure  above  the  plot.  The 
geometry  is  the  same  as  in  Figure  2.5  of  Chapter  2.  Units  of  RCS  are  not  used,  and  the 
frequency  is  300  MHz  instead  of  1  GHz,  so  the  levels  are  not  comparable  to  the  previous 
results  of  this  chapter.  Three  different  wind  speeds  are  considered  in  addition  to  the  flat 
surface  result  (wind  speed  zero).  For  5  m/s,  a  =  0.135  m  and  a/X  =  0.135.  For  10  m/s, 
a  =  0.54  m  and  a/X  =  0.54.  For  15  m/s,  a  =  1.22  m  and  a/X  =  1.22.  The  important  thing 
to  note  for  this  case  is  that  the  backscattered  field  can  peak  above  the  flat  surface  result  for 
all  three  wind  speeds. 
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Figure  5.3:  2D  RCS  vs.  time  for  a  13  m  high  ship  on  a  402.4  m  surface  at  three  different 
wind  speeds.  Elevation  angle  =  5°,  frequency  =  1  GHz. 


shipiand  surface 


Surface  #4  (wind  speed  =  15  m/s,  time  =  0  s  above) 


GFB  Tim 

e-evolvinc 

1  Surface  #4  - 

Coherent - 

Fiat  Surface  ■  • 

Time  (seconds) 


Figure  5.4:  2D  RCS  vs.  time  for  a  13  m  high  ship  on  a  402.4  m  surface.  Elevation  angle 
5°,  wind  speed  =  15  m/s,  frequency  =  1  GHz. 


Scat.  Field  Magnitude  (dB) 


30  .  Source  radiation  pattern  "'/■  30 
20  position 

10  ,y  £ 


MoM  Region 


0  20  40  60  80  100  120  140  160  180  200  x(m) 

Wind  Speed  =  15  m/s,  at  time  =•  0  s 
f=300MHz,  1=204.8,  angle=5  deg 


h . ». 


\;.s  ' 

t  j  \  ,/V> 


v/\ 

■:  t m v): 


■  I  W  t  1 

v  pr 

1  K:  / 


T  I  /, 


’  1  ' 
"  1 

U  •  :■ . I 


+  u=0m.s 

- u=5m/s 

—  u=10m/s 
- u=15m/s 


time  (sec) 


Figure  5.5:  2D  backscatter  vs.  time  for  a  13  m  high  ship  on  a  204.8  m  surface.  Elevation 
angle  =  5°,  frequency  =  300  MHz. 


5.3  RCS  as  a  Function  of  Ship  Roll  Angle 

Figure  5.6  shows  the  RCS  as  a  function  of  roll  angle  for  the  13  m  high  ship  on  a  flat  surface 
with  hull  angles  of  80°  and  90°  with  respect  to  the  horizontal  plane.  The  elevation  angle  of 
the  incident  plane  wave  is  5°.  The  figures  above  the  plot  show  the  ships  at  a  roll  angle  of 
+10°.  The  large  peak  around  10°  for  the  80°  hull  angle  is  due  to  the  left  hull  face  forming  a 
90°  corner  reflector  with  the  sea  surface.  Likewise,  the  large  peak  around  0°  for  the  90°  hull 
angle  is  due  to  the  left  hull  face  forming  a  corner  reflector  with  the  sea  surface.  The  peaks 
at  -5°,  5°,  and  15°  are  all  due  to  the  incident  or  sea  reflected  plane  wave  hitting  the  left  hull 
face  broadside.  Away  from  the  peaks,  the  pattern  variations  and  overall  levels  are  similar  to 
the  RCS  vs.  time  plots.  The  exception  is  the  90°  hull  angle  case  for  low  wind  speeds,  where 
the  corner  reflector  effect  dominates  independently  of  time. 

It  is  interesting  to  note  that  the  roll  angle  patterns  are  very  similar  for  roll  angles  of 
less  than  -5°.  For  these  angles  the  scattering  appears  to  be  dominated  by  the  ship  topsides 
(which  are  identical),  and  the  scattering  from  the  hull  has  a  lesser  effect. 
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RCS  (dBm) 


Flat  Surface,  Hull  Angles  =  80°  and  90°  (roll  angle  =  10°  above) 


Figure  5.6:  2D  RCS  vs.  roll  angle  for  a  13  m  high  ship  on  a  flat  surface.  Elevation  angle 
5°,  frequency  =  1  GHz. 
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5.4  Effect  of  Target  Geometry  on  the  RCS 


As  was  seen  in  the  last  two  sections,  the  angle  the  hull  makes  with  the  sea  surface  has  a 
strong  impact  on  the  ship  RCS  pattern.  A  hull  angle  of  90°  gave  the  highest  RCS  levels  for 
flat  and  time-varying  sea  surfaces.  In  general,  the  RCS  for  rough  time- varying  surfaces  was 
lower  than  the  flat  surface  result.  The  ship  geometries  considered  so  far  have  very  high  radar 
signatures.  Next,  the  scattering  characteristics  of  a  low-RCS  shape  will  be  investigated. 

Figure  5.7  shows  the  RCS  vs.  time  for  a  low-RCS  shape  with  wind  speeds  of  15  m/s 
and  5  m/s,  and  the  elevation  angle  is  again  5°.  Also  shown  is  the  flat  surface  RCS  for 
this  elevation  angle.  The  shape  is  9  m  high  and  the  hull  makes  an  angle  of  110°  with  the 
horizontal  plane  (20°  from  vertical).  The  overall  RCS  level  is  about  30  dB  lower  than  the 
previous  ship  geometries  (note  the  vertical  scale  is  shifted  with  respect  to  previous  plots), 
but  more  important  to  note  is  that  the  flat  surface  RCS  is  less  than  most  of  the  peaks  in  the 
time- varying  RCS  for  this  case.  It  is  also  significant  to  note  that  the  RCS  patterns  for  the  two 
different  wind  speeds  have  about  the  same  range  of  variation.  This  was  not  necessarily  the 
case  for  the  conventional  ship  geometries  considered  in  the  last  two  sections  of  this  chapter. 
In  particular,  the  RCS  vs.  time  plot  for  the  90°  hull  angle  at  a  wind  speed  of  5  m/s  was  very 
close  to  the  flat  surface  result,  while  the  15  m/s  result  had  large  variations. 

Figure  5.8  shows  the  RCS  as  a  function  of  roll  angle  for  the  same  low-RCS  shape  on  a 
flat  surface  (computed  using  coherent  scattering/image  theory).  The  elevation  angle  of  the 
incident  plane  wave  is  again  5°.  For  a  roll  angle  of  zero  the  RCS  has  a  minimum.  The  large 
peak  at  -15°  is  due  to  the  incident  plane  wave  hitting  the  left  hull  face  broadside.  At  -20° 
the  left  hull  face  is  vertical,  and  makes  a  90°  corner  reflector  with  the  ocean.  Since  under 
normal  conditions  it  is  unlikely  that  the  shape  would  roll  more  than  10°  from  vertical,  the 
expected  peak  RCS  level  is  about  10  dB  for  this  shape.  However,  the  RCS  vs.  time  plots  for 
this  shape  indicate  a  maximum  RCS  level  of  about  0  dB.  A  more  detailed  study  is  necessary 
to  determine  which  level  is  more  accurate  under  realistic  conditions,  but  the  conservative 
estimate  is  probably  the  most  advisable  to  follow  because  a  roll  angle  of  only  a  few  degrees 
greatly  increases  the  RCS  for  this  case. 
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It  is  clear  from  the  results  of  this  section  that  the  shape  of  the  target  has  a  significant 
effect  on  how  RCS  varies  as  a  function  of  the  time-changing  sea  surface.  It  is  clear  that  the 
roll  angle  may  be  the  deciding  factor  on  the  peak  RCS  level,  but  it  is  beyond  the  scope  of 
this  investigation  to  determine  how  the  roll  angle  of  a  ship  varies  with  the  sea  state. 

5.5  Conclusions 

The  RCS  vs.  time  results  of  this  chapter  have  shown  that  the  RCS  of  a  ship  on  a  rough 
surface  can  be  higher  than  the  RCS  on  a  flat  surface.  This  is  especially  true  for  low-RCS 
shapes.  On  the  other  hand,  the  RCS  as  a  function  of  roll  angle  for  a  ship  on  a  flat  surface 
shows  the  maximum  expected  RCS  values.  The  maximum  roll  angle  for  a  given  ship  and 
sea  state  should  be  available  in  most  situations  of  interest  (but  is  not  available  here) ,  so  the 
roll  angle  plot  is  probably  the  most  useful  for  determining  the  peak  RCS  level  which  may  be 
expected.  It  also  seems  to  provide  a  good  indication  of  the  expected  variation  of  the  RCS 
as  the  sea  surface  changes  with  time  or  with  wind  speed. 
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Chapter  6 

Conclusions  and  Future  Directions 

6.1  General  Observations 

The  2D  investigation  reported  here  into  the  radar  scattering  phenomenology  of  a  target  on 
a  rough  sea  surface  has  yielded  some  interesting  results,  some  expected,  some  not.  It  was 
expected  that  the  time- varying  RCS  of  a  ship  on  a  rough,  moving  sea  surface  could  have 
peaks  higher  than  the  RCS  of  the  same  ship  on  a  smooth  surface.  However,  this  appears  to 
be  true  mainly  for  low-RCS  shapes.  For  high-RCS  targets,  such  as  conventional  ship  designs 
where  the  hull  makes  close  to  a  90°  angle  with  the  sea  surface,  the  time-varying  RCS  is 
generally  lower  than  the  flat  surface  level.  For  low-RCS  targets  the  time- varying  RCS  ranges 
above  and  below  the  flat  surface  level. 

It  was  expected  that  the  variation  (i.e.,  maximum  deviation)  of  the  time- varying  RCS 
would  depend  strongly  on  the  wind  speed  (or  wave  height).  However,  somewhat  surprisingly, 
the  results  have  shown  that  the  variation  is  about  the  same  for  different  wind  speeds  (except 
for  the  case  when  the  wind  speed  is  very  low  and  the  target  has  a  high  RCS,  which  tends  to 
be  close  to  the  flat  surface  case). 

Although  a  great  deal  of  effort  was  spent  to  develop  the  spectrally  accelerated  GFB 
method  described  in  Chapter  2,  it  was  invaluable  for  efficiently  computing  many  reliable 
scattering  patterns.  This  made  possible  the  phenomenology  study  of  Chapter  5.  It  also 
provided  valuable  reference  solutions  for  evaluating  the  ray-optical  method  of  Chapter  4. 
Two  journal  papers  are  being  prepared  based  on  this  work. 
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The  investigation  of  Chapter  3  into  the  plane  wave  spectrum  of  the  sea  scattered  field 
(with  the  target  absent)  yielded,  as  expected,  a  singular  coherent  component  in  the  direction 
of  reflection  from  a  flat  surface.  In  addition  to  the  time-independent  coherent  component, 
the  incoherent  spectrum  was  seen  to  be  randomly  distributed  in  angle.  The  incoherent 
power  spectrum  appears  to  have  a  well-behaved  probability  distribution  in  angle,  whereas 
the  phase  is  uniformly  distributed.  Therefore,  the  incoherent  spectrum  averages  to  zero 
when  sampled  over  many  time  steps.  The  amplitude  of  the  coherent  component  of  the  sea 
scattered  field  is  well-predicted  by  the  physical  optics  formula  derived  in  [1].  It  approaches 
the  amplitude  of  the  incident  plane  wave  for  very  low  incidence  angles  and  wind  speeds,  but 
rapidly  approaches  zero  as  the  elevation  angle  or  wind  speed  increases,  or  the  frequency  of 
the  incident  field  increases.  A  discrete  plane  wave  expansion  was  also  investigated  Chapter  3, 
and  was  shown  to  adequately  reconstruct  the  incident  and  sea  scattered  field  in  the  vicinity 
of  the  target.  The  expansion  uses  the  incident  and  coherent  sea  scattered  plane  waves,  in 
addition  to  a  number  of  the  strongest  incoherent  components. 

The  plane  wave  expansion  for  the  incident  field  illuminating  the  target  was  used  in 
Chapter  4  as  a  ray-optical  excitation  for  finding  the  ship  scattered  field  over  the  surface  of 
the  ship.  As  an  approximation  which  allows  the  use  of  ray-optical  methods,  it  is  assumed  that 
the  ship  is  on  a  locally  flat  surface  when  solving  for  the  ship  scattered  fields.  A  reciprocity 
integral  was  derived  which  computes  the  target  backscattered  field  at  the  source  using  only 
the  incident  and  scattered  fields  over  the  surface  of  the  target,  without  having  to  back-track 
the  scattered  fields  back  to  the  source.  This  approach  effectively  decouples  the  ship  scattering 
problem  from  the  sea  surface.  It  was  found  that  by  using  only  the  coherent  components  of 
the  illuminating  field,  the  RCS  would  generally  be  under-predicted.  Using  a  more  complete 
discrete  plane  wave  expansion  using  50  plane  waves  predicted  well  the  level  of  the  RCS 
pattern,  but  not  the  details  because  of  the  locally  planar  approximation  for  the  sea  surface. 
However,  the  ray-optical  approach  developed  here  is  not  expected  to  predict  the  RCS  pattern 
exactly,  but  instead  be  able  to  predict  the  peak  levels  and  range  of  variation  of  the  pattern 
which  it  does  nicely. 
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The  discrete  plane  wave  expansion  used  for  the  incoherent  sea  scattered  field  was  obtained 
here  for  a  pre-determined  sea  surface,  and  required  a  numerically  intensive  computation.  It  is 
of  interest  to  investigate  the  possibility  of  using  a  randomly  generated  plane  wave  expansion 
based  on  the  probability  distribution  of  the  incoherent  spectrum.  This  will  be  part  of  future 
work. 

The  RCS  vs.  time  results  of  Chapter  5  were  very  insightful  for  visualizing  the  variation 
and  peak  levels  of  the  RCS  for  a  time-changing  sea  surface.  However,  since  these  results 
are  only  available  from  a  numerically  intensive  reference  solution  which  includes  the  entire 
sea  surface  (i.e.,  GFB),  it  is  desirable  to  be  able  to  instead  use  the  ray-optical  techniques 
which  remove  the  sea  surface  from  the  computation.  Unfortunately,  the  coherent  scattering 
computations  were  shown  to  generally  under-predict  the  RCS  level,  and  the  more  complete 
plane  wave  expansion  is  also  numerically  expensive  to  obtain  (although  this  may  be  remedied 
in  future  work). 

The  most  useful  ray-optical  results  were  shown  to  be  the  flat  surface  RCS  patterns,  both 
as  a  function  of  elevation  angle  and  roll  angle.  For  high-RCS  targets  the  flat  surface  result 
is  generally  the  worst  case  in  terms  of  peak  RCS  levels.  But  for  low-RCS  targets  the  flat 
surface  RCS  may  be  lower  than  the  peaks  of  the  time-varying  RCS.  In  those  cases  the  roll 
angle  plots  are  especially  useful  for  visualizing  the  variation  of  the  RCS  around  roll  angles 
close  to  zero.  The  flat  surface  RCS  is  easy  to  compute  using  ray  techniques  and/or  image 
theory. 

6.2  Future  Work 

Clearly,  there  is  a  little  more  investigation  needed  for  2D  problems  before  proceeding  to  3D. 
A  better  defined  plane  wave  expansion  is  needed  for  the  sea  scattered  field  which  illuminates 
the  ship,  one  which  only  needs  to  know  about  the  statistical  properties  of  the  sea  surface 
and  not  the  exact  surface.  The  expansion  should  use  the  fewest  plane  waves  necessary  to 
get  a  good  result,  because  too  many  plane  waves  could  greatly  reduce  the  efficiency  of  the 
ray-optical  method.  Another  possible  approach  which  is  currently  being  investigated  is  to 
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use  an  equivalent  planar  impedance  surface  or  periodic  surface  to  model  the  sea  surface 
locally  in  the  vicinity  of  the  ship.  The  equivalent  surface  would  be  chosen  to  best  simulate 
the  plane  wave  scattering  characteristics  of  a  random  sea  surface  as  a  function  of  wind  speed, 
frequency,  and  incidence  angle. 

The  extension  to  3D  will  begin  with  a  spectral  scattering  analysis  of  a  realistic  3D  ocean 
surface.  This  will  determine  the  statistical  properties  of  the  plane  wave  spectrum  so  that  a 
randomly  generated  discrete  plane  wave  expansion  can  be  used  to  illuminate  the  ship.  A  3D 
surface  scattering  code  based  on  the  forward-backward  method  with  spectral  acceleration  is 
currently  being  developed  under  a  separate  project. 

Next,  the  plane  wave  scattering  from  a  3D  target  will  be  computed  using  the  local  planar 
approximation  for  the  sea  surface.  The  scattering  will  be  computed  using  the  method  of 
moments  as  in  the  2D  case,  in  addition  to  high-frequency  methods  such  as  physical  optics 
and  UTD.  Of  particular  interest  will  be  the  use  of  the  RCS  Basic  Scattering  Code  (RCS-BSC) 
which  is  based  on  the  UTD  and  was  developed  here  at  the  OSU  ElectroScience  Laboratory 
[11]- 

It  will  also  be  necessary  to  obtain  3D  reference  solutions  for  a  target  on  a  rough  surface, 
possibly  again  using  numerical  solutions  for  some  limited  cases,  and  using  available  measure¬ 
ments.  Not  as  much  effort  will  be  put  into  a  full  3D  numerical  solution  (which  implicitly 
includes  the  sea  surface)  as  was  put  into  the  2D  GFB  code  because  of  the  numerical  and 
geometrical  complexity.  Only  a  few  selected  configurations  will  be  analyzed  numerically  to 
establish  the  reliability  of  the  ray-optical  methods. 

6.3  A  Final  Note  on  Coherent  Scattering 

It  needs  to  be  re-emphasized  that  for  most  realistic  situations  the  incidence  angle  is  very 
close  to  grazing  the  sea  surface  because  the  height  of  the  source  antenna  is  going  to  be  much 
smaller  than  the  distance  to  the  target.  Therefore,  the  coherent  sea  scattered  plane  wave 
component  is  expected  to  typically  be  very  strong  according  to  Eq.  (3.6).  Physically,  this 
means  that  the  random  variations  in  the  sea  scattered  field  tend  to  average  out  to  zero  over 
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long  propagation  distances,  so  only  the  coherent  field  remains.  This  lends  further  importance 
to  the  flat  surface  RCS  analysis  because  for  grazing  angles  the  incident  field  is  nearly  the 
same  as  for  a  flat  surface.  Unfortunately,  the  GFB  reference  solution  becomes  less  accurate 
for  grazing  incidence  because  the  sea  surface  must  be  truncated  to  a  finite  length  which 
is  usually  much  less  than  the  distance  from  the  source  to  the  target.  Since  the  numerical 
reference  solution  is  unreliable  for  this  case,  realistic  experimental  data  is  needed  to  support 
the  hypothesis  that  the  coherent  scattering  computation  is  sufficient  in  most  cases  of  interest, 
provided  that  variations  in  the  target  roll  angle  are  included  in  the  analysis. 

Figure  6.1  shows  RCS  vs.  time  plots  for  the  low-RCS  target  from  Figure  5.7  at  a  wind 
speed  of  5  m/s,  with  the  elevation  angle  of  the  incident  plane  wave  at  5°  and  1.6°.  The 
magnitude  of  the  coherent  sea  scattered  plane  wave  component  computed  from  (3.6)  for  an 
elevation  angle  of  1.6°  is  0.99  of  the  incident  plane  wave  magnitude.  Therefore,  the  RCS  is 
not  expected  to  vary  much  with  time  as  the  plot  confirms.  For  5°  the  coherent  component 
is  0.89  of  the  incident  plane  wave,  so  the  RCS  would  also  not  be  expected  to  vary  much 
with  time,  but  the  plot  shows  large  variations.  Investigation  of  realistic  experimental  data  is 
needed  to  better  understand  at  what  point  the  RCS  becomes  strongly  coherent,  and  whether 
the  time  variation  of  the  RCS  is  determined  more  by  the  sea  surface  or  the  rolling  of  the 
target.  One  thing  is  certain  from  the  numerical  results  of  this  report:  if  the  coherent  sea 
scattered  term  given  by  (3.6)  is  much  smaller  than  1,  the  time-varying  RCS  of  a  target  on 
the  sea  surface  will  have  large  variations.  A  limit  to  the  variation  is  reached  as  the  coherent 
component  approaches  zero,  which  depends  on  the  frequency  and  elevation  angle  as  well  as 
the  wind  speed. 
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